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PROPER 3Cs A NUCLEON-PION TRANSPORT CODE 


By Jules J, Lambiotte, Jr., John W. Wilson, 
and Tassos A. Filippas 
Langley Research Center 

SUMMARY 

A set of computer programs, referred to as the PROPER 3C Transport Code, is 
being used at the Langley Research Center to calculate and analyze the interaction of 
nucleons and pions within a prescribed medium. The programs simulate the transport 
of nucleons and pions (with energies up to 2000 MeV) within a slab. The necessary 
nuclear interaction cross sections are read from a master data tape which contains 
cross-section data for various elements. The nuclear interaction data above 400 MeV 
are incomplete at this time but will be updated as data become available. A particle 
biography tape is generated and used to compile distributions for particles involved in the 
transport. 

Descriptions of the programs, detailed user's guides, and sample input and output 
for programs in the code are presented. 

INTRODUCTION 


PROPER 3C is a nucleon-pion transport code which was developed at the Langley 
Research Center for use in radiation shielding studies. The bulk of this code is the 
PROPER 3B code given in reference 1. PROPER 3B was developed to simulate the 
transport of nucleons with energies below 400 MeV, PROPER 3C consists of modified 
versions of the programs in PROPER 3B, an additional program to extend the transport 
capability to 2000 MeV (2 GeV), and a replacement for the PROPER 3B statistical 
program. 

The purpose of this paper is to describe the capabilities of the code and to present 
a detailed user's guide. The theoretical considerations in the model are presented only 
generally since quite detailed explanations of the transport theory for the nucleons are 
available in the PROPER 3B description (ref. 1). The pion calculations, which are simi- 
lar, are described in the section entitled "PROBLEM DESCRIPTION." 

The computer programs in the PROPER 3C code use Monte Carlo techniques for the 
transport of neutrons, protons, and/or pions from energies as high as 2 GeV down to neu- 
tron thermal energy. Although the transport medium is restricted to having either a 



circular or square cross section, it may consist of a mixture of many elements. There 
are provisions to allow a source spectrum with any composition of the aforementioned 
particles, incident at any energy below 2 GeV and at any angle. A continuous slowing 
down model is used to slow the charged particles. Nuclear interactions also are allowed 
(presently only for nucleons since pion data are not available). A complete history of 
each source particle and its derivative particles is developed and recorded on a biography 
tape. This history includes the particle location within the medium, its direction of 
travel, energy, parent particle, and chronological generation. The information contained 
on this tape is used by a distributional analysis program to compute the distribution for 
each type of particle at any desired depth as a function of particle position or momentum. 

The organization of PROPER 3C is shown in figure 1. PROPER 3C contains three 
transport programs (TWOGEV, PROTOS, and SUPER B), a distributional analysis pro- 
gram (STAT), a master data tape which contains nuclear data for the transport, and two 
programs (INTPOL1 and INTPOL2) which rearrange the data from the master tape into a 
form usable by the transport programs. 


SYMBOLS 


d distance to nuclear interaction 

E^ neutron thermal energy 

f(x) probability density function for variable x 

Ij interval index for the jth variable in frequency distribution function (this func- 

tion is defined in description for program STAT) 

j ,K dummy variables 

LOG total number of computer locations reserved to store frequency distribution 

function 

LOCTOT number of computer locations needed to store frequency distribution function 
for all particles at any one depth 

N(ll >*2 >• •'’!)) frequency distribution function for j variables 

NB number of bits allotted to each storage bin for computing frequency distribu- 

tion function 
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NBIN 


number of bins to be packed into each computer word in calculating frequency 
distribution function 

NDEP total number of depths for which frequency distribution function may be cal- 
culated in one pass through biography tape(s) 

NLOC number of computer locations needed to store frequency distribution function 
for any one type of particle at any one depth 

NP number of different types of particles for which frequency distribution iunc- . 

tion is to be calculated 

n integer used to define multiplication in "splitting” or "Russian roulette" 

techniques 

nj number of intervals into which the range of the variable Vj has been divided 

R ,6r polar coordinates of projection of particle’s position onto X-Y plane (see 

fig. 2) 

Vj jth variable in frequency distribution function 

AVj interval length of variable Vj 

X,Y,Z Cartesian coordinate system for slab as defined in figure 2 

X',Y',Z' coordinate frame translated from X,Y,Z system such that particle is 
always positioned at the origin (see fig. 2) 

x generic variable denoting a member of a class of randomly distributed 

variables 

Xj[ ith member of a set of equiprobable values of x 

9,<p momentum angles defining direction of travel for a particle (see fig, 2) 

7r + ,7r°,7r“ pion triplet where + } o, and - refer to charge state 
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The following units are assumed throughout the paper unless otherwise noted; 


length cm 

energy MeV 

angle radians 


PROBLEM DESCRIPTION 

The problem to be solved is as follows; Given a source distribution of nucleons 
and pions (that is, their population, energy, and direction of motion) incident upon the 
face of a slab of material, determine the distribution of the particles within the slab. 

The problem is solved by using the Monte Carlo method. Each incident particle is 
transported separately through the specified volume. 

A continuous slowing down model, like the one described in reference 1, is used to 
slow the protons. In this model, the proton is slowed in discrete energy steps. The dis- 
tances the particle travels between these energy levels are obtained from a normal dis- 
tribution which describes the straggling. The mean and standard deviation of the normal 
distribution are computed from the theory of energy loss of a charged particle in matter. 
For charged pions, a continuous slowing down model also is used, but straggling is not 
computed. 

The path to nuclear interaction d is chosen from a set of random numbers 1 dis- 
tributed according to the exponential distribution function 

f(d) =^e“ d /“ 

where a is the collision mean free path. (See ref. 3, ch. 6.) When a nuclear collision 
occurs, the distribution of secondary particles is obtained from tabular data which have 
been entered into core storage from a master data tape. These data define the transport 
characteristics for each type of particle in the given slab material; that is, at various 
energy levels for an incident particle, the data include the mean free paths, the average 
number of particles produced from a collision during the cascade process and their 
energy and direction, and the average number of evaporation particles produced and their 
energy. The pions are not allowed to have nuclear interactions since data are not yet 

'‘"The random number generator used is the FUNCTION subprogram RANF supplied 
as software for the CBC series 6000 computer systems by Control Data Corporation. 

The characteristics of RANF are analyzed in reference 2. 
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available to define the secondary products. These data will be incorporated as soon as 
possible. For low-energy neutrons (below 18 MeY), data are entered into core which 
define the probabilities of absorption, neutron-neutron reactions, neutron-— two™ neutron 
reactions, fission, and elastic scattering. The data also define the angular distribution 
and energy after the reaction. The transport of the source particle is terminated if it 
leaves the slab, is absorbed, or suffers a collision. All particles resulting from the 
interaction are treated as secondaries and transported before the process is repeated 
for the next source particle. 

After the transport is completed, a final program, STAT, may be used to calculate 
desired distribution functions within the slab. 

The nuclear interaction data mentioned previously are stored in arrays at equiprob- 
able values of the independent variable (ref. 4). That is, if f(x) is the probability den- 
sity function for x, the K + 1 bounds for the K equiprobable intervals of x are the 
set xj (j = 1,2,...,K+1) where each x^ and satisfy the equation 



With this approach, a uniform random number is used to generate the interval from which 
x is obtained. 

Quite often when using the Monte Carlo method, one encounters difficulties because 
either too few particles survive at a certain depth to accurately describe the distribution 
there or because in certain regions many more particles survive than are needed for an 
accurate description. These two problems can be approached by using "importance 
sampling" techniques "splitting" and "Russian roulette,” respectively (ref. 3), Splitting 
means that in certain regions where insufficient populations of particles are present, one 
arbitrarily creates n particles for each initial particle present and assigns to each 
particle a weight of 1/n. The splitting technique should be employed only when the com- 
puter resources prohibit the use of more source particles and not even then if the colli- 
sion frequency distribution is not well defined. Russian roulette is the inverse of splitting 
and involves the transport of only one particle out of every n particles with a weight of 
n assigned to it. 

These two techniques are programed into the code. The user specifies the splitting 
and Russian roulette ratios. When a count of all the particles is made by program STAT, 
the weight assigned to each particle is used to determine whether the particles count as 
less than, equal to, or greater than a full unit. 



PROPER 3C DESCRIPTION 


Due to considerations of computer storage and execution time, the transport is 
conducted in three stages by the following three programs: 

(1) TWOGEV (Langley program R2092): TWOGEV transports neutrons, protons, 
and pions with energies up to 2 GeV, slows the protons to 400 MeV, slows the charged 
pions to any preselected cut-off, and considers nuclear interactions by the protons and 
neutrons. The pion nuclear interactions will be included when the data become available. 
Any nucleons created with energy below 400 MeV are stored on the output biography tape 
and treated, along with the protons which have been slowed to 400 MeV, by the next link 
in the code, PROTOS. 

(2) PROTOS (Langley program R2091): PROTOS treats protons below 400 MeV and 
neutrons in the 18 to 400 MeV range. These particles may be introduced as input data 
through a source spectrum or they may be obtained from a biography tape as a continua- 
tion of a TWOGEV run. The protons are slowed to a preselected energy and nuclear 
interactions are allowed. Any neutron created below 18 MeV is placed on the output 
biography tape and transported in the third link, SUPER B. 

(3) SUPER B (Langley program R2093): SUPER B transports only neutrons whose 
energy is less than 18 MeV. The neutrons are transported until they exit the slab, are 
absorbed, or lose energy down to a preselected cut-off. 

These three programs require nuclear data which include the kinds and probabili- 
ties of nuclear interactions and the number, energy, and directional distributions of 
derivative particles. This information is available for a limited number of elements 
(table I) on the master data tape. TWOGEV and PROTOS require auxiliary programs, 
INTPOL1 and INTPOL2, respectively, to prepare the data in final form. 

(4) INTPOL1 (Langley program R2090): INTPOL1 interpolates linearly in atomic 
weight between elements on the master data tape in the 2 GeV range and mixes the ele- 
ments to produce the desired shield. The data are output on a "Monte Carlo" tape for 
inppt to TWOGEV. 

(5) INTPOL2 (Langley program R2089): INTPOL2 performs the same function as 
INTPOL1 for data in the 400 MeV range. It generates a tape which is used in PROTOS. 

The particle distributions are compiled in the final link in the code, ST AT. 

(6) STAT (Langley program R2094): After the transport is completed, ST AT uses 
the particle-histories output on the biography tape to compute the particle distribution at 
points within the slab. 
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These six programs are described herein in the following order: First, the pro- 
gram INTPOL1 and INTPOL2, which supply the Monte Carlo tapes to the transport pro- 
grams; then, the three transport programs TWOGEV, PROTOS, and SUPER B; and, 
finally, the program STAT, The flow chart of PROPER 3 C is given in figure 1. 

PREPARATION OF MONTE CARLO TAPES 

Each of the three transport programs requires a nuclear data tape (Monte Carlo 
tape) as input to the program. The data to be used are on the master data tape which 
consists of three files containing the nuclear data for the elements specified in table I. 

The three files contain data in the energy range used by TWOGEV, PROTOS, and , 

SUPER B, respectively. For SUPER B, the data for a particular shield may be copied 
directly from the master tape if the data are present. (See table I.) However, TWOGEV 
and PROTOS require the auxiliary programs INTPOL1 and INTPOL2, respectively, to 
organize the data into usable forms. Abstracts and flow charts for these programs are 
provided in appendix A and appendix B, respectively. The data for TWOGEV are taken 
from Bertini's preliminary calculations (ref. 5) which assumed that pion production was 
through the (3, 3) -resonance and that pion decay was 25 percent forward, 25 percent back- 
ward, and 50 percent isotropic in the resonance rest frame. With the exception of the 
pion production, the model is the same as that described in reference 6. Comparison of 
Bertini's calculations with experimental data are given in reference 7. Reference 6 
describes the nuclear interaction model that generates the data used by PROTOS (that is, 
for energies below 400 MeV), and gives experimental comparisons. The low-energy neu- 
tron data for SUPER B has been taken from the United Kingdom Neutron Data File 
(ref. 1). 


Preparation of Data Tape for TWOGEV by INTPOL1 (R2090) 

Main program .- The main program is labeled INTPOL1 and defines the following 
tape (disk) storage units: 

TAPE9 — A disk storage unit onto which the user copies the appropriate file from the 
master data tape. It is later used for intermediate storage. 

TAPE? — A disk storage unit used for intermediate storage. 

TAPE 10 - A di.sk storage unit used for intermediate storage. 

TAPE8 - A disk storage unit which will contain the final data for all the shields at the 
completion of the run, 

TAPE5 and TAPES are the card input and printed output devices, respectively. 
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INTPOLl reads the nuclear data from the master data tape. It reads card input to 
describe the elements for which data are needed (see "Card input”). If data for the ele- 
ment specified are already available from the master data tape, INTPOLl merely copies 
the data to a disk file (TAPE 10). If the element data are not available, the program inter- 
polates linearly in atomic weight between the next lighter and heavier elements and puts 
these data on the disk (extrapolation is not permitted). As many elements as are needed 
to produce the desired shield or shields may be obtained in this fashion. Note from 
table I that only data for oxygen and aluminum are now available in this energy range so 
that the interpolation feature is of limited usefulness at this time. More elements will be 
added to the table as soon as data are available. 

Two additional tapes should be described at this point. TAPE9 is a disk file which 
contains the nuclear data for the elements noted in the TWOGEV portion of table I. 

TAPE 10 (disk) contains the interpolated data as an intermediate tape. The format of the 
two tapes which is identical is as follows: 


Record 1 


TAPE 9 

TAPE 10 

Z 

ZNEW 


Record 2 


TAPE9 

TAPE 10 

ATOM 

AA 

ELEMNT 

ELMENT 

BTAB 

BTABNY 

ATAB 

ATABNY 

NSPK 

NSPKNY 

NSPE 

NSPENY 

DSPK 

DSPKENY 

SIGMA 

S1GMNY 


The atomic number of the element 


The atomic weight 
Name of element requested 

Array containing the equiprobable energies of the 
evaporation particles in MeV 

Array containing the equiprobable cosines for the 
scattering angle 

An array containing the mean number of secondary 
cascade particles per reaction 

An array containing the mean number of secondary 
evaporation particles per . reaction 

An array containing the standard deviation from the 
mean of secondary cascade particles per reaction 

An array containing the total cross sections in barns 
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Records 3 through 7 

Each record is an array containing equiprobable 
energies of- a particular secondary cascade particle 
for each equiprobable cosine of the scattering angle. 
These data have been normalized by the incident 
energy value at which they were tabulated (see 
explanation of JE which follows) 

The arrays have the following effective dimensions: 

4 21 5 2 

BTAB (JE, I, IN, UT) 

4 11 5 5 

ATAB (JE, L, IN, UT) 

4 5 5 

NSPK (JE, IN, UT) 

4 5 2 

NSPE (JE, IN, UT) 

4 5 5 

DSPK (JE, IN, UT) 

4 5 

SIGMA (JE, IN) 

4 11 41 5 

FACT (JE, L, I, IN) 

where 



JE number of energy points for incident particle (the tables have been prepared 

at four energy values: 400 MeV, 750 MeV, 1000 MeV, and 2000 MeV) 

I number of values per IN-UT combination and energy group 

L number of cosine values per IN-UT combination and energy group 

IN type of incoming particle 

UT type of secondary particle 

The values of IN, UT are 

1 for proton 

2 for neutron 
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3 for 7r + 

4 for tt° 

5 for 7r“ 

It should be noted that the five records for FACT are data for UT = 1, 2, 3, 4, 5 in that 
order. These arrays are described in more detail in the BINOUT1 writeup, 

INTPOL1 then calls subroutine REDIG1 which combines the elements into the 
desired shield. Upon return from REDIG 1 one more revision is made to the data. Due 
to the large storage which would be required by a FACT array effectively dimensioned as 
FACT (4, 11, 41, 5, 5), the data must be packed. The data are packed six values to a 
word into a new integer array NFC which has been normalized so that the maximum value 
of the FACT array is stored in NFC with a value of 2*0 - 1 = 1023. The TWOGEV 
transport program unpacks the data when needed. 

The final data tape then has this form: 


Record 1 
IDENT 

Record 2 

DENSIT ^ 

ELMENT 

BTABNY 

ATABNY 

NSPKNY 

NSPENY 

DSPKENY 

LAMBLA 

PH 


See "Card input" (p. 13) 


As in the description in REDIG1 


Record 3 

C The multiplicative factor used to normalize the FACT array to 

1023 at its maximum 

NFC The FACT array packed six values to a word; NFC is dimensioned 

as NFC (7517) 

These three records are repeated for each shield. The data are stored on disk (TAPE8) 
as one file and should be copied to the user’s tape at the conclusion of the run by COPYBF 
control card. 
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Subroutines ,- INTPOL1 contains two major programs for combining data into proper 
combinations and format and for providing printed output, 

REDIG1: This subroutine uses the transport data from one or more of the elements 
on the intermediate disk (TAPE 10) to produce data for one or more shields which may 
have any specified composition (of up to five elements) or density. Data cards defining 
the composition of the shield(s) are read in. The transport data for each shield are gen- 
erated by using the data on TAPE 10 and are stacked one behind the other on the disk file 
designated as TAPE9. The data for each shield are in the exact format as the input tape 
to EED1G1 except that in the second record for each shield, the SXGMNY array has been 
replaced by two arrays 

LAMBLA mean free paths 

PH probabilities for hydrogen interaction 

These two arrays are dimensioned as 
4 5 

LAMBLA (JE, IN) 

4 5 

PH (JE, IN) 

where JE and IN have the same meanings as before. Data for hydrogen interactions are 
not yet included so that the user must not request data for compounds which contain 
hydrogen. A return is made to INTPOL1 where the data are packed and recorded on 
TAPE8 as previously described! 

BINOUT1: This subroutine, when requested, prints the contents of the tapes 
described previously. Three calls are automatically made by INTPQL1 to list the tapes. 
The first is to list the tape input to REDIG 1, the second is to list the tape output from 
REDIG1, and the third is to an entry point BIN1 in BINOUT1 which lists only the FACT 
array after it has been packed into and then unpacked from array NFC. This final print- 
out shows the effect of packing on the original FACT array data. Any of these listings 
can be suppressed by setting the appropriate sense switch 1, 2, or 3 by using the SWITCH 
n control card. The output is as follows: 

BTAB — This array defines the 20 equiprobable^ energy intervals for evaporation parti- 
cles. For every energy level there are two groups each having five columns and 21 rows. 
The five columns are for the five values of IN and the 21 rows are for the 20 equiprobable 
intervals. 

^The equiprobable values contained in BTAB, ATAB, and FACT have been com- 
puted by using the integral method mentioned in the PROBLEM DESCRIPTION. 
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ATAB — For each energy level there are 11 groups describing the 10 equiprobable cosine 
intervals for the cascade particles. Within each group there are five columns for the 
five types of particles in and five rows for the same five types out. 

NSPK - Within each energy level there are five columns for the five types of particles 
in and five rows for the same five types out. 

DSPK — This array has the same format as NSPK. 

NSPE - This array has the same format as NSPK except that there are only two rows 
corresponding to protons and neutrons out. 

SIGMA — There are four rows corresponding to the four energy levels and five columns 
for the five types of particles. 

FACT - This array defines for the cascade particles 40 equiprobable energy intervals 
determined at each of the 10 equiprobable cosine intervals (ATAB). The energy values 
are normalized to the incident energy levels. In this printed output there are five groups 
corresponding to the five types of secondary particles. Each group is divided into JE 
energy levels. Within each energy level there are five columns (corresponding to the 
five types of incoming particles) and 44 rows (11 values of equiprobable cosine L by four 
values of equiprobable final energy I). Due to the size of the FACT array, only the first 
four values of I are printed for each of the 11 values of L, The first four lines are for 
L = 1, the next four for L = 2, etc. 

For the printout of the REDIG1 output tape, LAMB LA (mean free paths), and PH (proba- 
bilities for hydrogen reaction) are on the tape instead of SIGMA (cross sections). They 
are printed under the same format as SIGMA. 

Tape input .- The first file from the master data tape should be copied to TAPE9. 
The master data tape should be unloaded before execution of the program. This can be 
done with the DROPFIL control card (see Sample Problem Control Cards for INTPOL1). 

Card input .- All input variables in this section and in the other input descriptions 
that follow are typed floating point or integer by the usual I, J, K, L, M, N FORTRAN 
integer convention unless indicated otherwise. 
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Card 

set 

Number 
of cards 

FORTRAN 

name 

Columns 

Format 

Description 



NUMBER 

1 to 5 

15 

Number of elements for which data are requested 

A 

1 







NSHILD 

6 to 10 

15 

Number of shields to be produced from the elements in the previous item 



ZO 

1 to 4 

A4 

Atomic number of the element in table I just lighter than the desired element 






(right adjusted within field) 


N 

AO 

5 to 14 

F10.0 

Atomic weight of the element in the previous item 


U 

Z1 

15 to 18 

A4 

Atomic number of the element in table I just heavier than the desired element 

B 

M 




(right adjusted within field) 


B 

V 

Al 

19 to 28 

F10.0 

Atomic weight of the element in the previous item 


R 

ZNEW 

29 to 32 

A4 

Atomic number of the new element Z0 = ZNEW ^ Z 1 * 



AA 

33 to 42 

F10.0 

Atomic weight of the new element A0 ~ AA < Al * 



ELEMENT 

43 to 46 

A4 

Alphanumeric name of the new element 



ZO 

1 to 4 

A4 

Atomic number of element 1 in shield 



NA 

5 to 14 

F10.0 

Number of nuclei per unit volume of element 1 multiplied by lO^ 4 



Z1 

15 to 18 

A4 

Corresponding data for element 2 (or blank) 



NB 

19 to 28 

F10.0 




Z2 

29 to 32 

A4 



o 

NC 

33 to 42 

F10.0 

Corresponding data for element 3 (or blank) 

: 

£t 

Z3 

43 to 46 

A4 



•JSJ 

ND 

47 to 56 

F10.0 

Corresponding data for element 4 (or blank) 

c 

JN 

s 

Z4 

57 to 60 

A4 

Corresponding data for element 5 (or blank) 


H 

NE 

61 to 70 

F10.0 


I 

I 

h 



Continuation of data for this shield on next card 

i 

D 

NH 

1 to 10 

F10.0 

Number of hydrogen nuclei per unit volume multiplied by lO^ 4 (NH must be 






set to 0.0 for INTPOL1 usage) 



DENSIT 

11 to 20 

F10.0 

Density of medium, cm3 (for identification only) 



NR 

21 to 25 

15 

Number of elements in the compound or mixture excluding hydrogen (e.g., for 






NaOH, NR = 2) 



IDENT 

26 to 29 

A4 

Identification number of the new material 



ELEMENT 

30 to 33 

A4 

Name or formula of the new material 


*If the new element ZNEW is in table I, the data should be merely copied to disk by setting ZO = Z1 = ZNEW 
and AA = AO ^ Al. 
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Tape output .- After execution, the final Monte Carlo data are contained in a file on 
intermediate disk (TAPE8). A blank tape should be mounted and the contents of TAPES 
should be copied onto it with a COPYBF control card. 

Printed output .- As the data for each element are read from TAPES in 1NTPOL1, 
the atomic number of the element is printed followed by the atomic number of the element 
requested (ZNEW). When REDIG1 has completed its function, it prints "REDIG OS," 

For the printed output of the contents of the tape, see subroutine BINOUT1 description. 

Preparation of Data Tape for PROTOS by INTPOL2 (R2089) 

Due to similarities between INTPOL1 and INTPOL2, only the differences between 
the two are noted in this section. The details can be found in the section describing 
INTPOL1. 

Main program .- The main program is labeled INTPOL2 and defines the same tape 
and disk units as INTPOL1. As seen in table I, the data for PROTOS (25 to 400 MeV 
energy range) is much more extensive so that the interpolation feature has more meaning. 
The formats of TAPE9 and TAPE10 which differ somewhat from those for INTPOL1 are 
as follows: 

Record 1 


TAPE 9 

TAPE 10 

Z 

ZNEW 


Record 2 


TAPE9 

TAPE 10 

ATOMV 

AA 

ELEMNT 

ELMENT 

BTAB 

BTABNY 

ATAB 

ATABNY 

FACT 

FACT 

NSPK 

NSPKNY 

NSPE 

NSPENY 

SIGMA 

SIGMNY 
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The effective dimensions of the arrays are 

9 21 2 2 

BTAB (JE, I, IN, UT) 

9 11 2 2 

ATAB (JE, L, IN, UT) 

9 11 41 2 2 

FACT (JE, L, I, IN, UT) 

9 2 2 

NSPK (JE, IN, UT) 

9 2 2 

NSPE (JE, IN, UT) 

9 2 

SIGMA (JE, IN) 

where JE, I, L, IN, and UT have the same meanings as before. The IN-UT indices are 
reduced in number because no pion data are used in PROTOS. The nine energy levels 
are for 25, 50, 100, 150, 200, 250, 300, 350, and 400 MeV. 

REDIG2 is called to combine the data for the shields. Upon the return to INTPOL2, 
FACT is packed as before into an NFC array but only four values to a word. NFC is 
dimensioned as NFC (4059). The final output tape (on TAPES, a disk file as before) is 
the same as the output tape from REDIG1 except that NFC is now a smaller array due to 
the fewer particle types to be considered. (Pions are not treated in PROTOS.) 

Again, three similar calls are made to BINOUT2 to provide printed output. 

Subroutines .- INTPOL2 contains two major programs for combining data into proper 
combinations and format and for providing printed output, 

REDIG2: The function of REDIG2 is identical to that of REDIG 1. The output tape 
has the same format as the input tape except that SIGMNY has been replaced by 
LAMBLA (9, 2) and PH (9, 2) as before. 

Cross sections for hydrogen interactions are contained in a data statement so as to 
include shields containing hydrogen. 

BINOUT2: Although similar in function to BINOUT1, the output from BINOUT2 is 
simplified because there are only four IN-UT combinations and these may be put into four 
labeled columns, 

ATAB — Within each energy level there are 11 rows defining the 10 equiprobable cosine 
intervals. 
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BTAB — Within each energy level there are 21 rows defining the 20 equiprobable energy 
levels. 

FACT — Within each energy level there are 44 rows whose meaning is described under 
BINOUT1. 

NSPK — There are nine rows corresponding to the nine energy levels. 

NS PE — This array has the same format as NSPK. 

SIGMA — There are two columns for neutron and proton cross sections and nine rows for 
the nine energy levels. 

LAMBLA — This array has the same format as SIGMA. 

PH — This array has the same format as SIGMA. 

Tape input ." The second file from the master data type should be copied to TAPE9. 

Card input . - It is the same as the input of INTPOL1 except that NH may be nonzero 
since interactions with hydrogen are possible. 

Tape output . - It is identical to the output of INTPOL1. 

Printed output .- It is identical to the output of INTPOL1. 

Nuclear Data for SUPER B 

The nuclear data for SUPER B need only to be transferred from the third file of the 
master data tape to TAPE9 by means of control cards. This file contains the low-energy- 
neutron data for the elements listed in table I. There are two records per element. The 
COPYBR control card is used to put the two records for the desired element onto the 
disk TAPE9. The user is restricted to these elements because there is no interpolation 
capability. 


NUCLEON-PION TRANSPORT 

The transport is carried out through three programs: TWOGEV, PHOTOS, and 
SUPER B. 


TWOGEV (R2092) 

Main program . - The main program is labeled TWOGEV and defines the following 
tape (or disk) storage units: 

TAPE 4 — The biography tape. 

TAPES — The source tape (or disk). St is either input to the program or generated 
within it. 
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TAPE9 — The Monte Carlo data tape. 

TAPE5 and TAPES are the card input and printed output devices, respectively. 

The purpose of TWOGEV is to call various subroutines and to terminate the output 
biography tape at the conclusion of the run. 

Subroutines .- The main program, TWOGEY, performs mainly the managerial task 
of calling subroutines which read input, prepare data, perform the calculations, and ter- 
minate the output. 

PDATA: This subroutine reads from cards the stopping power data used in the 
program to simulate the slowing down of the protons and computes the tables necessary 
to define this process. (See ref. 1.) Data from the nuclear data tapes are read into core 
also. 

PIDATA: As an entry point under PDATA, it is similar in function to PDATA 
except that the input pertains to the pions. Although the charged pion rapidly decays into 
a muon, no distinction is made between the two in the slowing down process since the 
stopping power data for the two is very similar (the difference being about 5 percent). 

GDATA: This subroutine reads from cards data defining the geometry of the 
transport medium. No calculations occur in this subroutine. 

GEO: The transport is executed in this subroutine. A particle is taken from the 
source tape that was generated in the source subroutine. The type of particle determines 
the action to be taken. A proton may either have a nuclear interaction or merely be 
slowed in energy. A neutron has only nuclear interactions. The charged pions lose 
energy through continuous slowing down. They are presently not allowed to have nuclear 
interactions because of the lack of data. Each charged pion decays weakly in 10"® sec- 
ond (half-lifetime) into a muon and neutrino (approximately 99 percent of the time). The 
neutrino is ignored and the muon is slowed down. The neutral pion tt ° decays electro- 
magnetically in 10" *6 second into two gamma rays. Because of this short lifetime, the 
7 r°'s are only counted and recorded on the biography tape. If the particle is to be slowed, 
the distance it will travel before losing energy to the next energy level is first taken from 
the slowing down tables and then the particle is transported that distance and its energy 
is reduced to the next level. The particle is slowed continually until it either has a 
nuclear interaction or until it reaches a preselected cut-off energy. For a particle which 
is to have a nuclear interaction, its free path to interaction is chosen from a set of ran- 
dom numbers distributed according to the exponential collision frequency distribution, 
and the particle is moved to its interaction site. A call is made to subroutine GENU 
which uses the nuclear data to determine the kinds and number of secondary particles 
and places this information in a buffer. GEO then takes another particle either from this 
buffer or, when it is empty, from the source tape. 
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For a proton which may either interact or slow down, its free path and slowing 
down distance are compared, with the lesser distance determining the action taken. This 
test is made at each energy level in the slowing down tables. 

GENU: When a particle suffers a collision, subroutine GENU uses the necessary 
nuclear data to determine the number and kinds of secondary cascade and evaporation 
particles, their energies, and their directions. 

TEMP: This subroutine is called from GENU to facilitate the conservation of 
energy in the nuclear collisions (actually only the fourth component of the four-momentum 
is conserved). A temporary storage is set up to hold data for all the secondary particles 
resulting from one nuclear interaction, TEMP has an entry point TRY. 

TRY: After all secondary particles have been accounted for, the sum of their 
energies is compared with the energy of the original particle minus the binding energy 
loss. If the comparison is not sufficiently close, a return is made to GENU and the pro- 
cess is repeated. When a good comparison is found, the pertinent data for the secondary 
particles are stored in a buffer in the BUFOUT subroutine. 

BUFOUT: This subroutine contains a storage buffer where data pertaining to the 
secondary particles can be stored and retrieved. 

The data are stored by a call from TEMP to an entry point BUFIN, The data 
stored for each secondary particle consists of the particle type, position, direction of 
movement, energy, energy index within the slowing down tables, relative weight, parent 
particle, and generation. The buffer is initially set through a call to ENTRY BUFSET. 

A maximum of 100 particles may be stored at any one time. Any attempt to overstore 
results in a program exit. 

The data are retrieved from the buffer by a call to BUFOUT from GEO. The last 
particle put in the buffer is retrieved first. Should the buffer be empty, a flag is set to 
inform GEO to take the next source particle instead. 

SOURCE: This subroutine defines the source spectrum incident on the face of the 
shield. Although any source spectrum made up of the types of particles considered in 
the transport may be included by the user, certain standard spectra are already pro- 
gramed into the subroutine and may be requested. These include neutron and/or proton 
source spectra that are normally or isotropically incident, are either spread uniformly 
over the surface or all incident at the origin, and have energies uniformly distributed 
over a specified interval. 

The source spectrum is written on logical unit 8 (TAPES) by the following call to 
subroutine PUT which is a tape input-output routine: 

CALL PUT (8, MS, X, Y, Z, SLAT, CTH, STH, CFI, SFI, NP, E, NSOU) 
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where 

MS index of the energy S in the slowing down table SK; that is, MS is smallest 

integer such that EK(MS) = E (see description of EX that follows) 

X,Y,Z particle position (see fig. 2) 

SLAT relative weighting of the particle normally set to 1 for primaries 
CTH,STH cosine and sine of the angle 6 (see fig. 2) 

CFI,SFI cosine and sine of the angle <j> (see fig. 2) 

NP particle type (1, 2, 3, 4, 5 for proton, neutron, ir + , ir°, it", respectively) 

E energy of the particle 

NSOU source particle identification number 

When the user writes his own source program, he has the following variables 
available, most of which are read into the computer in subroutines PDATA and GDATA: 

EK(100,2) energy values for the slowing down tables; the second index is 1 for proton 
data and 2 for pion data 

M(2) number of values in the proton and pion tables (a maximum of 100) 

NT number of source particles 

RU radius or half-side of the transport medium 

DU length of medium, 0 < Z < DU 

IG =1 for square-cross-section transport medium 

= 0 for circular-cross-section transport medium 

EO highest point on energy scale and usually the source energy for a mono- 

energetic spectrum 

These variables are mentioned again under ’’Card input." 

PUT: This input-output subroutine uses the RECIN and RECOUT library subrou- 
tines and calls DPACK and UNPACK. 

DPACK: This assembly language subroutine packs two computer words into one 
location. 

UNPACK: This assembly language subroutine unpacks one computer word into two 
locations. 

ST ATI: This subroutine is called to determine whether the particle is still within 
the limits of the transport medium. It also stores the particle data on the biography tape. 
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A record is written on the biography tape (logical unit 4, TAPE4) by a call to sub- 


routine PUT: 


CALL PUT 

(4, X, Y, Z, SL1, SL2 

,, EPR, IK1, NSOU, SLAT, IGEN) 

where 



X,Y,Z 

particle coordinates 


SL1,SL2 

angles 6 and 0 


EPR 

particle energy 


IK1 

integer defining the particle type and status, which may have the following 
values: 


I, 2, 3, 4 
-1,-2, -3, -4 

II, 12,13,14 
-11,-12,-13,-14 

0 

10000 

|lKl| modulo 10 

for protons 
for neutrons 
for ir+ 
for 7T“ 
for v° 

to denote the end of a case — that is, to denote that all 
particles in the source spectrum have been 
transported 



= 1 for a moved particle which is still inside the 
boundary 

= 2 for a moved particle which has passed outside the 
boundary 

= 3 for a particle moved to a collision site 
= 4 for a new particle 

NSOU 

parent particle identification number 

SLAT 

relative weighting of the particle 

IGEN 

chronological generation of the particle 


Three things are done to reduce the number of biography tapes. First, PUT calls 
a COMPASS subroutine DPACK which packs two values into one word so that there are 
only five words per record. Second, RECIN and RECOUT are used to block the output. 
Third, the program is allowed to write to the end of the tape and, thus, there is no wasted 
space. A DAYFJLE message will note the EOT (end of tape) condition and the next tape 
loaded will also be entered. No data are lost with this procedure. ST ATI has several 
entry points. 
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OUT; As an entry point under ST ATI, it is called to output particle data on the biog- 
raphy tape, 

OUTL; As an entry point under STATI, it is called when the run terminates to write 
the final record on the biography tape. 

CHAN1L: This subroutine computes the elements necessary to transform from a 
system that moves with the particle to the stationary X,Y,Z system as described in 
reference 1, 


Tape input .- TWOGEV expects the Monte Carlo tape generated from INTPOL1 to 
be mounted defined as TAPE9. A blank tape should be mounted defined as TAPE4 to be 
the output biography tape. 

Card input .- The main program and several subroutines require card inputs. 

The main program TWOGEV controls the flow of the program by means of an input 
variable I which indicates the subroutine to call. The variable I is read with an 15 
FORMAT and may have the following values: 

1=1 calls PDATA, returns; then calls PEDATA and returns. 

1 = 2 calls GDATA and returns. 

1 = 3 calls GEO (which calls SOURCE) and returns after transport completed. 

1 = 6 reads 80 characters from the card immediately following this one. 

This heading is output from GEO when that subroutine is called. 

I = 7, 8, 9, or 10 terminates the run; an end of file is written on the biography tape and a 
program exit is made. 

The subroutine PDATA contains a NAMELIST, PD AT, which includes the following 
variables: 

EO largest energy to be included in the proton slowing down tables (should be 

equal to or greater than the energies of the source particles) 

ES cut-off energy for protons 

N number of points in the stopping power tables 


EN 

DES 

M 


N energy values where stopping power tables are generated, ENj < ENj + i 


corresponding N values of stopping power, 


MeV 

g/cm^ 


number of energy points desired in proton slowing down tables 
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KE indicator for the kind of energy scale in the slowing down tables 

KE = 1 gives M - 1 logarithmic steps between EO and ES 

KE = 2 gives M - 1 linear steps 

KE = 3 indicates that the slowing down table is read instead of computed 

EK table of slowing down values; not input unless KE = 3 

RO density of shield, gy/cm3 

ZE effective atomic number of shield 

CF straggling parameter taken from figure 3 

UNIF a starting value for a uniform distribution of random numbers 

RAN a starting value for a normal distribution of random numbers 

AV atomic or molecular weight 

The subroutine PIDATA contains a NAMELIST, PIDAT, which includes the following 
variables: 

EP largest energy to be included in the pion slowing down tables 

EPS cut-off energy for pions 

NP number of points in the pion stopping power tables 

MP number of energy points desired in the slowing down tables 

EN same as in PDATA but for pion data 

DES same as in PDATA but for pion data 

The subroutine SOURCE is called from GEO to prepare the source tape for the 
transport. A NAMELIST SOURCE contains the following variables: 

ISO = 1 for all source particles incident at (X = 0, Y = 0, Z = 0) (see fig. 2) 

= 2 for source particles spread uniformly over the Z = 0 plane 
= 4 for source program to be written by the user 

This indicator generates a transfer to the end of this subroutine where 
the user may insert his FORTRAN IY cards. 

IS I =1 for normal incidence 

= 2 for isotropic incidence 

FAC fraction of neutrons in source spectrum 

El lowest energy in source spectrum 

E2 highest energy in source spectrum (for monoenergetic source El = E2) 
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The subroutine GDATA contains a NAMELIST, GDAT, which has the following 


variables; 


IG 

= 1 for square -cross-section transport medium 
= 0 for circular-cross-section transport medium 

DU 

length of transport medium 

RU 

half-square side or radius 

NT 

number of source particles 


Tape output .- The particle histories are recorded on the biography tape defined as 

TAPE4. 

Printed output .- The main program TWOGEV prints the indicator I. 

Subroutine PDATA prints the input data and the slowing down table. 

Subroutine PIDATA prints the corresponding pion data. 

Subroutine GDATA prints the variables input under GDAT. 

Subroutine GEO immediately prints the header card. After the transport for this 
case is completed, the results are printed in seven columns. The first column indicates 
the chronological generation of the particles for that row. Columns 3 to 7 contain the 
number of protons, neutrons, w+, v°, and tt~, respectively, in the generation indicated 
in column 1. The value in column 2 of each row refers to the data in columns 3 to 7 of 
the next row. It is the number of particles from the generation which suffered nuclear 
collisions. 

Subroutine SOURCE prints the source energy interval and the composition of the 
source spectrum. 


PHOTOS (R2091) 

PHOTOS and TWOGEV are nearly the same in both composition and usage. There 
are three main differences. First, PHOTOS transports only nucleons with energies less 
than 400 MeV ; whereas TWOGEV transports both nucleons and pions with energies up to 
2000 MeV. Second, PROTOS may use a biography tape from a TWOGEV run as a source 
tape; whereas TWOGEV must generate its own source tape. Third, to improve the sta- 
tistical results, a splitting procedure is incorporated in PROTOS whereby the user may 
have all source particles split at their first nuclear interaction. That is, for an n to 1 
splitting, when each source particle suffers its first nuclear interaction, GENU is called 
n times instead of once and each of the derivative particles carries a weight of 1/n 
instead of 1. This weight must then be accounted for in STAT. No such procedure is 
programed into TWOGEV. 
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The following description of PROTOS contains only new information since many of 
the subroutines are identical to their TWOGEV counterparts. 

Main program . - The main program is labeled PROTOS. The program card defines 
the tapes INPUT, OUTPUT, TAPE 4, and TAPE 9 as in TWOGEV. There is an additional 
unit, TAPE?, on which may be loaded a biography tape generated in TWOGEY to be used 
as a source tape. 

Subroutines .- As in TWOGEY, the main program PROTOS calls subroutines which 
read input, prepare data, perform the calculations, and terminate the output. 

PDATA: There is no entry PIDATA. 

GDATA: See TWOGEV description. 

GEO: This subroutine transports only neutrons and protons. For n to 1 splitting, 
it calls GENU n times. Under directions from the user, it takes source particles from 
a previously generated biography tape, TAPE7, instead of the usual source tape, TAPE8. 

GENU: See TWOGEV description. 

TEMP: Total kinetic energy is conserved in this version of the program. 

BUFOUT: The maximum number of allowable particles to be stored is 50. 

SOURCE: This subroutine is the same as in TWOGEV except that NP = 1 or 2, EK 
is dimensioned as EK (100), and M is not dimensioned at all. 

PUT: See TWOGEV description. 

ST ATI: The variable denoting particle status IK is restricted to neutron and proton 
values of IK1 only. 

CHANIL: See TWOGEV description. 

DPACK: See TWOGEV description. 

UNPACK: See TWOGEV description. 

Tape input .- PROTOS expects the Monte Carlo tape generated from INTPOL2 to be 
defined as TAPE9. If PROTOS is to be a continuation of a TWOGEV run, the output biog- 
raphy tape from that run should be defined as TAPE?. A blank tape should be mounted as 
TAPE4 to be the output biography tape. 

Card input . - The PROTOS card input is the same as the TWOGEV card input except 

that 

(1) There is no entry point PIDATA and no namelist P3DDAT. 

(2) ISO may be 3, in which case the source particles are taken from a previously 
generated biography tape defined as TAPE7, 
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(3) Namelist GDAT contains an additional variable WHT which is the inverse of the 
splitting desired. For example, to split 10 to 1, WHT = 0.1. Each particle with all its 
derivatives is assigned a weight (SLAT) equal to WHT. (See STAT for suggestions in 
using this option.) Variable WHT is set to 1.0 internally before GDAT is read. 

Tape output .- The histories of the particles are recorded on the tape defined as 

TAPE4. 

Printed output .- The PROTOS output is similar to the TWOGEV output except for 
obvious differences resulting from the fewer types of particles considered. 

SUPER B (R2093) 

Very few modifications have been made to SUPER B at Langley. For that reason, 
most of the following description consists of excerpts from reference 1. 

The PROTOS program is not useful at lower neutron energies where more detailed 
data are needed. If PROTOS creates a neutron at a lower energy than 18 MeV, it does 
not try to continue the transport process but puts it on the biography tape. The SUPER B 
program reads and copies the biography tape and when it finds a neutron with energy 
below 18 MeV, it transports the neutron and its derivative neutrons until they are either 
absorbed, leave the slab, or have a collision. Any neutron with energy below a prese- 
lected cut-off value is ignored. For any neutron with energy below thermal energy (a 
value specified by the user), SUPER B uses a constant mean free path and absorption 
probability. The program puts the biography of each event onto the new biography tape 
so that the final result is an expanded tape in the PROTOS format. The transport of the 
neutron is done by the standard Monte Carlo technique, that is, the mean free path at this 
energy is found, the actual path length is selected from the exponential distribution, and 
the particle is moved to the new point. If the neutron is still inside the transport medium, 
the absorption probability is compared with a random number to test whether absorption 
occurs. If absorption does not occur, the probabilities for different reactions are found 
and the element and reaction are selected. The resulting energy and scattering angle are 
determined and transformed to the laboratory system. The polar angle is transformed to 
direction cosines and the entire process is repeated for the particle at the new energy or 
for a new particle. 

Main program ,- The main program is labeled SUPER B and defines the following 
tape (or disk) storage units: 

TAPE 4 — The input biography tape. 

TAPE3 - The output biography tape. 

TAPE 10 — A disk storage unit onto which the user has copied the two appropriate records 
from the third file of the master data tape. 
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TAPE5 and TAPE6 are the card input and printed output devices, respectively. 

SUPER B reads the nuclear data into core and calls subroutine SUPB. 

Subroutines .- The main program SUPER B calls subroutines which read input, pre- 
pare data, perform the calculations, and terminate the output. 

SUPB: Subroutine SUPB reads the card data input, sets up initial values, and calls 
subroutine START. 

START: Subroutine START reads the PROTOS biography tape until it finds a neu- 
tron of energy less than 18 MeV. Since it sometimes happens that in a problem one gen- 
erates far more neutrons in the 0 to 18 MeV range than are needed for statistical pur- 
poses, an option is included which allows the user to employ a Russian roulette procedure 
to reduce the number of neutrons transported and hence the computation time and the 
number of output tapes. The Russian roulette procedure is the inverse of the splitting 
procedure described in PROTOS. In SUPER B, the user specifies the fraction of low- 
energy neutrons he wishes transported, with each one having its weight divided by this 
fraction. The neutrons which are not transported are merely deleted from the biography 
tape. By using the specified fraction and a uniform distribution of random numbers, it is 
determined whether this particle is to be transported. If so, computation proceeds to 
ENTRY EXPEN. 

EXPEN: This subroutine computes mean free path, computes the new particle 
position, and tests the absorption probability. If the neutron survives, the control goes 
to REAKT. If the particle is absorbed, control goes to START. EXPEN uses STAT1 
and OUT to put information on the biography tape. 

REAKT: This subroutine uses the probabilities entered as part of the nuclear data 
to choose first the atomic element of the shield and then the type of reaction. If the 
reaction is fission, the energy of the new particle and the laboratory scattering angle are 
computed in REAKT and control goes to REOR. If the reaction is not fission, control is 
given to the appropriate reaction of ELAST (elastic scattering), INELAS (n,n'), or 
NTVN (n,2n). 

ELAST: This subroutine tests the neutron energy to decide whether isotropic scat- 
tering occurs in the center-of-mass (c.m.) system. If not, it selects an angle from a 
given distribution. It converts the c.m. angle to the laboratory system and calculates the 
resulting energy in the laboratory system. Control then goes to REOR. 

INELAS: This subroutine tests the neutron energy to see if continuous or discrete 
energy loss is to be applied. If discrete, a level is chosen from the probabilities and the 
resulting laboratory energy is computed. The scattering is isotropic in the c.m. system. 
If continuous, the correct parameters for this element are chosen, an energy is drawn 
from the distribution, and control goes to REOR. 
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NTVN: TMs subroutine gets its secondary energy with an even chance from one of 
two distributions of the same general form as in the continuous part of INS LAST. 

RE OR; This subroutine tests the new energy to find out whether it is below cut-off. 
If so, control goes to START to search for a new particle. If not, REOR calls CHANIL 
and when control returns to RSOR, it calls EXPEN. 

CHANIL ^ 

STATI 

PUT > See PROTOS description. 

DPACK 

UNPACK 

Tape input . - The main program SUPER B expects the two records of interest from 
the third file of the master data tape to be copied to the disk unit designated TAPE 10. 

The biography tape containing the data for the low-energy neutrons from either a 
TWOGEV or PROTOS run should be mounted as TAPE3. The user must request a blank 
tape defined as TAPE 4 to be the output biography tape. 

Card input .- All input variables are defined through a namelist INPU as follows: 

EMIN cut-off energy 

DU length of transport medium 

RU radius or half-side of square 

IG =1 for square-cross-section transport medium 

= 0 for circular-cross-section transport medium 

NRECPT the first NRECPT records which go on the biography tape for this case will 
be printed 

IRAN random number 

Ell incident energy interval for this case 

E2 J (for identification purposes only) 

ETH thermal energy 

PAAS probability of absorption at ETH 

ALTH mean free path at ETH 

WHT fraction of low-energy neutrons to be transported, 0 < WHT = 1 

Tape output .- The histories of the particles are recorded on the biography tape 
defined as TAPE4. 
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Printed output .- The input data are printed and the NRECPT records are printed 
in columns as X, Y, Z, Q, <j>, E, IK, NSOU, SLAT, IGEN. 

The total number of records on the output tape is printed, 

COMPILATION OF PARTICLE DISTRIBUTIONS BY STAT (R2094) 

STAT is the final link in the PROPER 3C code. From the histories previously 
recorded on the biography tape, STAT obtains the information required to compile fre- 
quency distribution functions for particles in the transport. Frequency distributions are 
computed for a particle at any depth within the transport body as a function of the parti- 
cle's Cartesian or polar coordinates, energy, and/or direction of travel. The output may 
be in tabular or graphic form and/or on magnetic tape. 

Frequency Distribution Function 

The frequency distribution function N(li,l2,...,Ij) which is computed in STAT may 
be best explained by first looking at the one -dimensional case. Consider the flow of a 
particular type of particle across a plane at a depth Z. Let the range of a variable Vi 
be divided into nj equal intervals AVj. Then, the frequency distribution function 
N(li) for this type of particle at a depth Z within the body is defined as the number of 
times the variable Vi, when observed for the particles at this depth, was found to have 
a value lying in the Ijth interval; that is, 

(il - l) • AVi < Vi < Ii • AVi 

The extension to a function of j variables N(l^,l2,...,Ij) having nj,n2,...,nj 
intervals of lengths AVi,AV2,...,AVj is obvious. The values of N(li,l2,...,Ij) will be 
effectively stored as a j -dimensional array, and for future reference an array element 
will be called a bin. A maximum of five variables may be included in the distributions. 
These variables are particle coordinates X and Y (or R and energy E, and 

momentum angles 6 and <£. 

The user must specify the geometry of the volume used for the transport, the depths 
(values of Z coordinates) at which he desires the counting to be done, and the types of 
particles he wishes counted. STAT computes N(li,l2,...,Ij) by counting in one of two 
ways. The first method counts the number of particles crossing a specified plane 
(Z = Constant). If the path of a requested particle intersects one of the planes, the values 
of the variables are obtained at that depth by linear interpolation and the appropriate bin 
is updated. The second method involves counting through a volume. The user specifies 
another plane at Z + AZ which defines a subslab of the original body or slab, STAT 
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uses the specifications on X and Y (or R and GrJ to further partition the subslab 
into smaller volumes. Now each time the path of a particle intersects one of the bound- 
aries of these volumes, the appropriate bin is updated by 1/2. Figure 4 illustrates a 
typical situation. The particle intersects bin 4 at A, the face of the subslab, intersects 
the plane common to bin 4 and bin 2 at B, and goes out the back of bin 2 at C. Bin 4 gets 
updated by 1/2 for A, both bins 2 and 4 get updated 1/2 for B, and bin 2 gets updated 1/2 
for C. 

If one or both of the coordinates are not requested as variables, their interval is 
effectively taken by the program to be to +°° (o to 360° for For instance, if Y 

were not requested, there would be no intersection at B in figure 4. 

Calculation of Storage Requirements 

A one-dimensional storage area LC is set aside in STAT to hold the frequency 
distribution function and update it as the biography records are read from 

tape. The dimension of this array can and should be modified according to the needs of 
the particular computer run to make most efficient use of the computer. To be able to 
specify this dimension, the user should be acquainted with the method of storage and com- 
putation used in STAT. First, the storage area is used to compute N(li,l2,...,Ij) as a 
function of all the variables which appear in any of the desired distributions. Then, to 
compute the particular distribution N(l Ik), where K = j, which is a function of a 
subset of the variables in the previous section, a summation is performed over the 
unneeded variables. For example, suppose N(X,Y) and N(X,0) are requested; 

N(X,Y,0) is computed and later summed as follows: 

N(X,Y) =^N(X,Y,0) 

e 

N(X,6) =^N(X,Y,0) 

Y 

The storage area is constructed as follows. The number of locations NLOC required 
to store the function for any one type of particle at any one depth is calculated from the 
user's input specification. For NP types of particles requested, the number of loca- 
tions needed per depth is 

LOCTOT = NP*NLOC 

Now if LOC locations are set aside for the storage area by the user, the total 
number of depths for which N(lj,l2,...,Ij) may be calculated in one pass through the 
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biography tape(s) is determined from NDEP = LOC/LOCTOT by using integer arith- 
metic. In each pass through the biography tape(s), N(li,l 2 ,...,Ij) is calculated for 
NDEP more depths until all the depths have been processed. Obviously, the minimum 
amount of residence time, but also the maximum program length, occurs when enough 
locations are set aside to accommodate all the depths desired in just one pass through 
the tapes. The trade-off between the two must be made by the user. The modification 
to the program can be done easily (1) by changing the dimension of the LC array in the 
COMMON statement of the main program and (2) by setting LOC equal to this dimension 
as the first executable statement of the program. 

In order to decrease the storage required, the frequency distribution function 
N(ll,l2, •••>!]) is packed NBIN bins to a computer word. The user specifies NBIN and 
must be careful in the selection to avoid overflow of a bin. The number of bits NB 
which are allotted to each bin is computed from NB = 59/NBIN by using integer arith- 
metic, and thus the maximum value a bin could have would be 2^ B - 1. A message is 
printed which tells the number of times an overflow was attempted. A bin is not allowed 
to overflow. Its value remains at the maximum. The value of NLOC, needed previ- 
ously, may be calculated from 



NBIN 

by using integer arithmetic. The LC array is used as follows: Starting at LC(1), the 
first NLOC locations contain N(lj,l2,...,Ij) for particle type 1 at depth 1, the second 

NLOC locations contain the corresponding information for particle type 2 at depth 1, 
etc., up to a maximum of the four types of particles that can be requested. The same 
format for depth 2 is followed starting at LC(LOCTOT + 1). This format continues for 
as many depths as are being processed in this run. 

Thus, from the given values for the indices Ii,l 2 ,...,Ij and by knowing NLOC 
and LOCTOT, ST AT can compute the LC element LOCI for the NGth type particle 
at the Nlth depth by (1) computing LOC, the position of this particular element for the 
first particle type at the first depth, and then by (2) modifying this number as follows: 

LOCI = LOC + (NI - l)*LOCTOT + (NJ - l)*NLOC 

Then LC(LOCl) would contain the desired bin. 
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Weighting 


As stated in the PROTOS and SUPER B descriptions, the particles may be weighted 
differently due to splitting in PHOTOS and Russian roulette in SUPER B. This difference 
in weighting is accounted for in STAT by an integer input NORM, Since the packing 
scheme used requires an integer count in storage, each bin must be updated by an inte- 
gral amount. The variable NORM should be the smallest integer value which when mul- 
tiplied by each of the possible weightings will make all the products integers. Before 
output of results, each unpacked bin is divided by NORM to yield the correct distribution. 
Since NORM has the effect of reducing the maximum count a bin can have by a factor of 
NORM, one can make both good and bad choices for the weightings in the transport. 

As one example, the user splits 7 to 1 in PROTOS and keeps 0.3 of low-energy neu- 
trons in SUPER B. The possible weighting factors from PROTOS are 1 and 1/7 and from 
SUPER B are 1/0.3 and (l/7)/0.3. Therefore, all possible weights are 1, 1/7, 10/3, and 
10/21 and NORM must have a value of 21. 

As a second example, the user splits 7 to 1 in PROTOS and keeps 0.25 of low- 
energy neutrons in SUPER B. The possible weighting factors from PROTOS are 1 and 
1/7 and from SUPER B are 1/0.25 and (l/7)/0.25. Therefore, all possible weights are 
1, 1/7, 4, and 4/7 and NORM must have a value of 7. 

Clearly, the choices in the second example are better than those of the first exam- 
ple by a factor of 3 even though the weightings are nearly the same. 

Main Program 

The main program is labeled STAT and defines the following tape (or disk) storage 

units: 

TAPE4 — The input biography tape. 

TAPE 10 — An output tape which will contain the frequency distribution functions if tape 
output is requested. 

TAPE5 and TAPE6 are the card input and printed output devices, respectively. 

All card input is read in the main program. The logic necessary to identify the 
variables in the statistics is formulated. The storage requirements for the frequency 
distribution function are computed and subroutine GEOM is called to calculate the func- 
tion for the first case and for as many depths as allowed by the storage. Subroutine 
SUMM is called to compute the desired statistics by summing over the unwanted vari- 
ables. GEOM is called repeatedly until all the cases have been processed. If there are 
more depths, the tape is rewound and all cases processed again for these depths. 
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Subroutines 


Program STAT calls the following subroutines in order to read the biography tape 
and to compute and output the frequency distribution function for the variables specified 
by the user in the main program. 

GEOM: This subroutine, with the aid of subroutines INTSEC, COUNT, ORDER, and 
LOCAT, obtains the frequency distribution function for one case. It reads the biography 
tape until it finds a type of particle which has been requested and then traces the move- 
ment of the particle through the body by reading successive records from the tape. The 
location of the depths or subslabs to be processed has been supplied from the main pro- 
gram. If the path of the particle intersects any of these planes or subvolumes (fig. 4), 
INTSEC is called to obtain the data at the intersection point. Subroutine COUNT is then 
called to determine, when counting through a volume, whether there are other points of 
intersection. When IK = 10 000 is encountered on the tape, GEOM returns control to 
STAT. 

INTSEC: Given particle data at two points from the biography tape, subroutine 
INTSEC computes by linear interpolation the value of the variables at some intermediate 
point. 

COUNT: Once GEOM has determined that the path of the particle intersects any 
part of the plane or subslab, COUNT is called. If this run is to count across the plane, 
ORDER is called and given the variables at the intersection point from which it deter- 
mines the indices Ij,l2,...,Ij described previously. When counting through a volume, 
COUNT tracks the particle path by computing the intersections within the subslab (fig. 4) 
and supplies these values to ORDER. 

ORDER: Subroutine ORDER uses the values of each of the variables given from 
COUNT and the user's input defining the variables to compute the indices I|,l2,...,Ij. It 
calls LOCAT to determine the appropriate bin and then UPDAT to update that bin. 

LOCAT: Subroutine LOCAT uses the index values supplied by ORDER to compute 
which LC array element contains the bin and in which part of the word the bin lies. 
LOCAT has entry points RTREVE, which unpacks and retrieves this part of the word, and 
UPDAT, which updates the bin by some integral amount. 

INIT: Subroutine INIT computes values used by LOCAT in unpacking and updating 
the LC array. 

SUMM: An ordering array is constructed in SUMM to relate the position of each 
variable in the particular distribution to its position in the larger distribution which is a 
function of all the variables. The particular distribution is then computed by summing 
over the unneeded variables in the larger distribution. To retrieve N from the LC 
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array for a particular set of indices, SUMM calls AD which makes use of the ordering 
array. The retrieved values are adjusted for the weighting factor and output in a tabular, 
graphic, or tape format. After this distribution is computed and output, a return is made 
to STAT. 

PREPLT: Subroutine PREPLT substitutes Display Code characters for the values 
in the two-dimensional number density function to be output as a plot. 

PLOTN: Subroutine PLOTN plots N^I^. The range of N is divided into 50 equal 
intervals and the user's specification for the variable defines the abscissa. The symbols 
P, N, +, and - are used to denote for protons, neutrons, n + , and v~, 

respectively. 

AD: Subroutine AD uses the ordering array from SUMM to supply LOCAT with the 
correct ordering of the indices. The value of the bin is retrieved from RTREVE and 
return to SUMM. 


Tape Input 

The biography tape should be defined TAPE4. If there is more than one tape, the 
labeling provided by the operators at the time of the tape generation should handle any 
continuation problems. 


Card Input 

The card input for STAT is as follows: 
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Card 

set 

Number 
of cards 

FORTRAN 

name 

Columns 

Format 

Description 



NPLNS 

1 to 5 

15 

Number of depths at which N is to be computed 



NVAR 

6 to 10 

15 

Total number of different variables which will appear in any of the requested distributions (NVAR ^ 5) 



NSTAT 

11 to 15 

15 

Number of different distributions required (NSTAT = 10) 

A 

1 

IOPT 

16 to 20 

15 

= 1 to count across a plane 
= 2 to count through a volume 



NCASES 

21 to 25 

15 

Number of cases on biography tape(s) for which this set of input is applicable (NCASES = 20) 
(Each case terminates with a record having IK = 10 000.) 



NTAP 

26 to 30 

15 

Number of reels containing the particle data 



TAP 

41 to 50 

A10 

Tape number of first reel, left justified (not needed for NTAP = 1) 



NPOP(l) 

1 

11 

— 1 if want distributions for neutrons 
= 0 if do not want distributions for neutrons 



NPOP(2) 

2 

11 

Same as in previous item for protons 

B 

1 

NPOP(3) 

3 

11 

Same as in previous item for PI- 



NPOP(4) 

4 

11 

Same as in previous item for PI+ 



NBIN 

5 to 10 

16 

Number of bins to be packed into one word 



NORM 

11 to 15 

15 

Smallest normalization factor which makes all possible relative weights perfect integers 



XM(1) 

1 to 10 

F10.2 

Maximum value of X (or R) 

C 

1 

XM(2) 

11 to 20 

F10. 2 

Maximum value of Y (or blank) 



XM(3) 

21 to 30 

F10.2 

Maximum value of Z (see fig. 2) 


A 

Z<1) 

1 to 10 

F10.2 

First depth at which distribution desired 


s 

DZ(1) 

11 to 20 

F10.2 

Depth of first subslab (set to 0.0 if IOPT = 1) 


m 

a 

n 

Z(2) 

21 to 30 

F10. 2 

Second depth 

D 

a 

s 

DZ(4) 

71 to 80 

F10. 2 

Depth of fourth subslab 


n 

e 

e 

d 

Z(5) 

1 to 10 

F10.2 

Fifth depth 


e 

d 

DZ (NPLNS) 


F10. 2 

Depth of NPLNS subslab 



VN 

1 to 5 

A5 

Variable name (left justified) 

(See text that follows for possible entries.) 

E 

N 

V 

A 

D 

NBLK 

6 to 10 

15 

Number of intervals for this variable 

(Note restrictions on % and <f> in text that follows.) 


It 

V 

11 to 20 

F10. 2 

Minimum value of first bin* 



DELT 

21 to 30 

F10.2 

Change in value for each bin* 



NU 

1 to 5 

15 

Number of variables in this distribution 



VST(l) 

6 to 10 

A5 

Name of first variable in the distribution (left justified) 
(See text that follows for possible entries.) 



VST (2) 

11 to 15 

A5 

Name of second variable in the distribution (or blank) 


N 

VST(3) 

16 to 20 

A5 

Name of third variable 

F 

S 

T 

VST (4) 

21 to 25 

A5 

Name of fourth variable 


A 

T 

VST (5) 

26 to 30 

A5 

Name of fifth variable 



TAIP 

31 to 35 

15 

= 0 for no tape output of this distribution 
= 1 for tape output 



PLOT 

36 to 40 

15 

= 0 for no plot of this distribution 
= 1 for a plot 

No plotting can be done for NU > 2 

G 

N 

C 

A 

El 

1 to 10 

F10.4 

Lower bound of incident energy spectrum 


S 

E 

S 

E2 

11 to 20 

F10.4 

Upper bound of incident energy spectrum 


When variable is an angle, then input unit is degrees. 
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The names of the seven variables which might appear in card E or F (in ¥N or VST) 
are X,Y,R, THER, THEP, PHIP, or E which are for the coordinates X and Y (or R 
and momentum angles 6 and <p, and energy E. Due to the rotational symmetries 

of the variables % and cp about the Z-axis and the line Y = X in the X-Y plane, 
their values are always reduced by ST AT to a range of 0 to 45° and the bins should be 
allocated accordingly. 

This entire set of data cards is assumed to be pertinent to the NCASES run. If 
more distributions are sought, an entire new set of cards is required. The biography 
tape is always in a rewound position when card A is used. 

Tape Output 

The frequency distribution function N(li,l2,...,Ij) may be recorded on magnetic 
tape if requested. The tape should be mounted on logical unit 10. The data will appear 
on tape in the same order as they appear in the printout. For each case and at each 
depth, the following information is recorded per distribution: 

Record 1 

Z,DZ(NPOP(I),I = 1,4),E1,E2 
where these variables have been previously defined 

Record 2 

NUV,(LP(I),WS(I),DELS(I),LF(I),I = 1,NUV),NR 

where 

NUV the number of variables in this distribution 

LP(J) an array element whose value indicates the variable to which the Jth element 
of the other arrays refer 

LP(J) = 1 for X 

= 2 for Y 

= 4 e 

= 5 (f> 

= 6 E 

= 7 R 

= 8 % 


VVS(J) the lower limit of the first bin in the variable indicated by LP(J) 
DELS(J) the change in value from bin to bin for this variable 
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LF(J) the number of bins allotted to this variable 

NR an integer which tells how many records follow to give N for one type of 

particle 

Records 3 through 3 + NR - 1 


These NR records contain the frequency distribution function N for the first type 
of particle requested. Each record contains the following information: 

NWS , (INDX(I) ,1 = 1 ,NW) , (SUM (I) ,1 = 1,NWS) 


where 

NWS the number of elements in the SUM array 

NW NUV-1 

INDX indices of first NW variables 

SUM an array whose Ith element gives N for fixed indices INDX and the Ith 

interval in the last variable in the statistic 


There follows NR more records for each type of particle requested. 

The same type format is used on the remaining distributions. Next, the whole pro- 
cedure is repeated for the next depth and the following depths. Then the procedure is 
repeated for the next case on the biography tape. 

An example of a typical tape output is as follows: 

Presume that N(0,E) is requested for protons and neutrons at Z = 10.0 and AZ = 0.0. 
The incident energy was from 150 to 200 MeV, The momentum angle 6 has two bins in 
range 0 to 90° and the energy E has two bins in range 0 to 200 MeV. 

Presume also that the following information was compiled by ST AT: 

Neutrons Protons 


e 

E 

0 to 100 

100 to 200 


0 to 45 

18 

5 

45 to 90 

15 

8 


e 

E 

0 to 100 

100 to 200 


0 to 45 

25 

13 

45 to 90 

16 

11 


Then the output tape would contain 
Rec 1 10., 0., 1, 1, 0, 0, 150., 200. 

Rec 2 2, 4, 0., 45., 2, 8, 0., 100., 2, 2 
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Rec 3 

2, 1, 

25, 13 

Rec 4 

2,2, 

16, 11 

Rec 5 

2, 1, 

lO 

CO 

’S—i 

Rec 6 

2,2, 

15, 8 


Printed Output 

Information which describes the input is printed out. The maximum value which 
any one bin may have is printed as well as the number of attempts made to overflow the 
bins. After reading all the data from one case, printed out are IK, the total number of 
records, and the number of low neutrons. If IK ^ 10 000, there has been a read error. 

The distribution function N is recorded in one or more of the following forms 
(one always gets the tabular form except as noted): 

Tabular 

The range and interval size for each variable in the particular distribution is 
printed out. N(li,l2,...,Ij) is printed out by listing the indices of the first j - 1 vari- 
ables and printing N as the jth index varies over its range. 

Graphs 

Graphs may be obtained for NU < 3. For NU = 1, a plot of N^Ijj as a function 
of I is given on the printed output by using the PLOTN subroutine. For NU = 2, a 
rectangular grid based on the variable intervals is established and N is printed out by 


II 

O 

Blank 

N = 1 to 9 

1 to 9 

N = 10 to 35 

A to Z 

N > 35 

* 


The values for the asterisks are printed out below the graph. For this plot the 
tabular output is suppressed. 

Restriction : If a plot for N(li,l2) is requested, there is no tabular printout nor can 

there be tape output if I 2 > 120. 

SAMPLE PROBLEM 


A sample problem consisting of a source spectrum of 500 particles (40 percent pro- 
tons and 60 percent neutrons), isotropically incident on the top of the atmosphere with 
energy uniformly spread over the interval 1500 SE = 2000 MeV, was solved by means of 
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PROPER 3C, The transport medium is an infinitely wide slab of air 1000 cm thick. The 
air density was assumed to be 1 g/cm3. Because of the lack of high-energy data for 
nitrogen, data for pure oxygen have been used for air above 400 MeV. The stopping power 
data were obtained from reference 8, A splitting of 2 to 1 is performed in PROTOS and 
50 percent of the low-energy neutrons are transported in SUPER B. A total of approxi- 
mately 20 min of central processing time and 60 min of peripheral processing time was 
required to complete the transport. 

Listings of control cards, data input cards, and printer outputs are given for the 
sample problem. 
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SAMPLE PROBLEM CONTROL CARDS 



SAMPLE PRIBLfc M CONTROL C4RDS FOR INTPQLi 


CARD REFERENCE 


J0B 9 U300*7100. 

RUMCSI 

REQUEST TAPE1UHY* 20601 7 * RQ e JJL ^MASTER DATA TAPE 

RFW I NO ! T A PEI 1 I 

COP YRF! T APF 11 9 T APE9 ) 

DRDPFILCTAPEIH 
REWIND! TAPE9) 

SWITCH 3* 


L SO* 

REQUEST DIJMM2 * HY . 212054*RIL*JJL®MC 2GEV AIR 

REWIND! TAPE 8 ) 

COPYBF(T4PE8*DUMM2) 

UNLOAD! DUMM21 
EXIT® 

UNLOAD! DU MM2 I 
7/8/9 CARD 


JOB CARD 
COMPILE 

MOUNTS MASTER TAPE 

PUT HIGH ENERGY DATA 
ON DISK 

UNLOAD MASTER TAPE 

SUPPRESSES BINOUT1 
PRINTOUT OF PACKED 
FACT ARRAY 

MONTE CARLO TAPE 
FOR TWOGEV 

COPY DATA TO MC TAPE 


SAMPLE PROBLEM CONTROL CARDS FOR INTPOL2 


CARD REFERENCE 


JOB, 1*300*74000. 

RUN! S> 

REQUEST T APE1 1 * HY« 206017* RO,JJL ,M\ STER DATA TAPE 
RFWlND(TAPEil) 

COPYBF(TAPEIUDIS) 

C3PYBF!TAPEU*TAPE9) 

DR3PFI L { T APE1 1) 

REWIND! TAPE9) 

LSD. 

REQUEST 0UMM2 ® HY. 1202019 * Rl L* JJL*MC DATA AIR 

REWIND! TAPES! 

COP YBF ! TAPES* DUMM2 ) 

UNLOAD! DU MM2) 

EXIT® 

UNL 3AD ( DUMM2 ) 

7/8/9 CARO 


MOUNTS MASTER TAPE 

SKIP FILE l 
COPY FILE 2 


MONTE CARLO TAPE 
FOR PROTOS 

COPY DATA TO MC TAPE 
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SAMPLE PROBLEM CONTROL CARDS FOR TWOGEV 


CARO REFERENCE 


JDB*i # 400* 72DOO* 

RUN ! S I 

REQUEST DUMM2*HY» 2 1 2054, RQ • JJL * MC 2GEV AIR 
C3PYBF (DUMM2*TAPE9# l ) 

REMIND! TAPE9) 
rwOPFIL ( DUMM? } 

RFQUFST T APE4®HY * 50901 1 * R I L • J J L 9 T 4 OGEV TEST CASE 2 
LGO. 

UNLOAD I TAPE4I 
FX I T a 

UNLOAD! TAPE4J 
7/8/9 CARO 


SAMPLE PROBLEM CONTROL CARDS FOR PROTOS 


JOB*1*2000*63000<, 

RUN I S I 

REQUEST DUMM? *HY • 1202019 , RU . JJ L ,MC DATA AIR 
COPYBF(OUMM2*TAPE9*n 

REW I NO! TAPE9) 

DRQPFtL { DUMM2 ) 

REQUEST T APE7 » NY. 50901L *R0* JJL, TWOGEV TEST CASE 
REQUEST T APE4 $HY ® SAVTP* RIL»JJL«P*OT TEST CASE2 
LGO. 

UNLOAD! TAPE7I 
UNLOAD! TAPE4) 

FXI T® 

UNLOAD! TAPE7 ) 

UNLOAD! TAPE4) 

7/8/9 CARD 


MONTE CARLO TAPE 
PLACE HIGH ENERGY 
AIR DATA ON OISK 


OUTPUT BIO TAPE 


CARD REFERENCE 


MONTE CARLO TAPE 
PUT 400 MEV AIR DATA 
ON DISK 


SOURCE TAPE 
OUTPUT BIO TAPE 



SAMPLE PROBLEM COMTROL CAROS FOR SUPER B 


jrm*i*60Q»47oo* 

RUNtSS 

REQUEST TAPPlhHY. 206017 , RO* JJL * MA STER DATA TAPE 
C7PYBF(TAPPll*DI S s 2 ) 

C7PYBR4TAPPU,DIS»24> 

C7PYBRI TAPPll »T APE 10*2 ) 


DROPFIL(TAPPU) 

REQUEST TAPE3* HY® SAVTP* RIL 9 JJL*SB TEST CASE2 
REQUEST TAPE4* HY* 211050* RU*JJL, PROT TEST CASE2 
LSOo 

UNLOAD? TAPE3) 

UNLOAD I f A PEA-) 

FX I T* 

UNLOAD! TAPE3 I 
UNLOAD? TAPE4) 

7/8/9 CARD 


SAMPLE PROBLEM CONTROL CAROS FOR STAT 


J 7 B * l * 400 » 6 50 00 ® 

RUNIS1 

REQUEST T APE4*H Y* 1 07021 * RO * JJL * SB TEST CASE2 
REQUEST T APE1 0# NY® S AVTP , R IL * J J L * STAT TEST CASE2 


LGO* 

UNLOAD! TAPE4* 
UNLOAD I TAPEIO ) 

EXIT* 

UNLOAD ITAPE4) 
UNLOAD ITAP610) 
7/8/9 CARD 


CARD REFERENCE 


MOUNTS MASTER TAPE 
SKIP TWO FILES 
SKIP 24 RECORDS 
PUT LOW ENERGY 
NEUTRON AIR DATA 
ON DISK 

OUTPUT BIO TAPE 
INPUT BIO TAPE 


CARD REFERENCE 


INPUT BIO TAPE 
REQUESTED ONLY IF 
WANT TAPE OUTPUT 
FROM STAT 
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SAMPLE PROBLEM INPUT DATA 
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SAMPLE PROBLEM CARO INPUT FOR INTPOLX 


CARD REFERENCE 


l 

8 

i 

16*0 

8 

18 * 

8 L 6* 0 

AIR 

CARD 

CARD 

SET 

SET 

A 

B 

8 

0* 

*03764 

i* 


1 

8 « AIR 


CARD 

SET 

C 



SAMPLE 

PROBLEM 

CARD INPUT FOR INTP0L2 

CARD 

REFERENCE 

2 

1 



CARD 

SET 

A 

6 

12*0 

8 

16* 7 14* N2 

CARD 

SET 

B 

8 

16*0 

8 

18® 8 16® 02 




7 * 

0344 

8 * 

00752 

CARD 

SET 

C 

0* 

l* 


2 7*2 AIR 





44 



SAMPLE PROBLEM CARD INPUT FUR TWQGEV 


CARD REFERENCE 


6 

ATM3SPHERI C CALCULATION 503 PARTICLES ISOTROPIC 

1 

SPOAT EO=2000.0»R(l=l,000»N=74* M = 5 3 ? ES-400 ^0 ZE - 2 e C F=-5 * ? 
UMIF=®3,RAN=.3,AV=14*4*EN( 1 ) =2 . 0* 4*0#6*0, 8*0, 10.0 , 14 *0, 1 8.0, 22.0 * 26, 0* 
30.0, 34.0,38*0, 42.0,46*0, 50*0*60.0, 70. 0#80®0, 90 *0,1 00.0, 110.0,1 20.0* 
130.0* 140.0, 150.0,160.0* 170.0, 183® 0* 190 .0 *200*0® 220* 0*240*0* 260*0, 

2 80® 0 » 300 ® 0® 320.0*340.0# 360.0# 380. 0 *400*0 ,420. 0*440* 0,460 *0*480® 0* 

500*0, 520 * 0,540*0, 560 ,0« 580*0* 603®0s620®0,640*0*660®0*680®0 9 700*0 # 

720*0, 740*0, 760*0,780*0,800*0,8 23* 0,840*0,860*0*880*0,900*0,920*0® 
940*0,960*0, 980*0# 1000*0, 1200® 0 , 14 00* 0* 1600* 0# 2000® 0 , 

OESIl) =141 *7 #84 *16# 60® 22, 

4®727E*-0i,3.9920E*01,3®049E*0l,2o491E-*-Ql,2®118E-fr0t,i.851E$-01,l®650EQl, 
U492E01* 1.364E01, l*259EOl,1.17iEOl*1.097E01,9.5QO,8.428,7.608,6.959# 
6*432, 5# 997* 5® 629 * 5® 3 16 * 5® 054* 4® 839*4*601 ,4.417,4. 2 52*4. 104* 3® 970* 

3*738* 3® 54 3* 3. 378, 3 . 237, 3. 113, 3 .00 6,2 . 91 l ,2. 826,2 . 75 1 , 2 . 683 ,2. 622 , 

2* 567, 2* 516* 2® 471, 2.428,2* 393, 2. 35 4, 2* 322*2,291* 2® 263,2.237, 2.213, 

2® 190*2.169*2.149*2*130*2® 113* 2.097*2® 081 #2*067,2® 053*2,041*2*029, 
2*017,2.006, 1 * 996# 1*987# 1*978, 1.96 9*1 *961 * 1 » 899, 1.861, 1*839# 1® 819$ 

%? I DAT EP=2000» 0* NP-=74#MP = 50,EPS=2* 0* 

DES (11=31.50, 18. 02, 13. 01, 10. 35, 8. 697, 6* 730, 5® 597, 4. 8 57, 4. 336# 3. 949, 
3.651, 3*414, 3.222,3*062,2*929*2.673,2*492,2.358,2*256*2® 177,2.114,2.062, 
2.021*1.986, 1.958*1*934, 1.913, 1.896,1.882*1*870,1.851,1.838,1.828,1.822, 
1.818,1.816, 1*816,1.816*1.818, 1*320,1.823,1.826,1.830*1*834*1.838, 1.842* 
1.847, 1. 852, 1*857,1*861# l. 866# 1*871,1.876,1.881,1.886,1*891*1*896,1*901, 
1.906* 1.911, 1.916*1.920, 1.925, 1.930,1.934,1.939,1.944,1.948*1.953,1.957, 
1.999,2.037,2.071,2*130$ 

2 

$GDAT DU=iO0O*,RU=5OO0. , I G= 1, \)T = 500$ 

3 

&SORCE ISO=l,lSl=2,FAC=.4,Ei=15O0. #£2=2000. $ 

7 


1=6 FOR HEADING 
HEADING 

1=1 TO CALL POATA 
NAMELIST PDAT 


NAMELIST PIDAT 


1=2 TO CALL GDATA 
NAMELIST GOAT 
1=3 TO CALL GEO 
NAMELIST SORCE 
1=7 TO TERMINATE RUN 
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SAMPLE PROBLEM CARD INPUT FOR PROfOS 


CARD REFERENCE 


6 

ATMOSPHERIC CALCULATION TEST CASE 
1 

&PDAT EO=400.00®RQ=l,000*N = 80*M = 50*ES=l,000*KE=2*ZE = 7.,2*CF=-5.® 

UN i F =. 3* RAN=. 3 *AV=14* 4 ®ENU >=2.0*4. 0*6, 0*8.0* 10.0® 14.0*18.0*22.0*26.0* 
30.0® 3 4* 0® 38«0»42®0»46#0?50«0® 5 0» 0 f70*0?80»0j 90 ©0*100,0*1 1 0 © Q * 1 2 0 © 0 * 
130.0*140.0* 150*0® 160*0*170*0* 180, 0* 190.0*200. 0® 220. 0*240*0* 260.0® 

2 80*0® 300,0® 320,0*340,0® 360, 0® 380* 0 *400.0*420. 0* 440. 0*460 .0*480.0 ® 
500,0® 520, Oe 540,0*560,0*580,0*500* 0*620,0*640.0*660,0*680,0*700,0* 

7 20, 0» 740, 0* 760, 0® 780 *0 * 800, 0 ® 8 20. 0*840 ,0 * 860.0* 880 * 0* 900 ,0 *920,0® 
940.0*960,0*980*0* 1000,0® 1200, 0 * 1 4 00, 0 * 1600, 0 * 2000. 0*2400 ,0 *2800. 0 * 
3200.0* 3600, 0*4000,0® 50OO.O®DE5(l)=l*4i7Es-O2g>8*416E-s-OX®6.O22E*Ol» 

4. 727E*0l® 3* 9920E4- 01 * 3. 049EF0 1*2. 49 lE^-01 *2.11864-01® l . 85 IE *01 * i . 650EO l * 
U492E01*1.364E0i® 1.259E01® 1.171601*1,097601*9.500*8,428*7.608*6,959* 
6,432* 5.99 7*5.629*5.316® 5, 054® 4.80 9,4. 60 1*4*417 *4* 252® 4. 104*3,970® 

3. 7 38 »3» 54 3*3,378*3.2 37 *3. 113*3. 00 6*2. 911*2. 826*2. 751*2,6 83 *2.622® 

2» 56 7®2»516»2. 471*2.428*2,390,2, 354*2.322*2.291*2.26 3* 2. 2 37 *2. 213* 
2.190*2* 169*2. 149*2.130* 2,113*2.097*2,081 *2.067*2.053*2.041*2,029* 
2.017*2.006® 1.996® 1.987® 1.978* 1 .96 9 * 1 , 961 * 1 . 899® 1 . 86 1 * 1 . 8 39® 1 * 819 * 
1.816* 1. 822® 1.833® 1.846® 1.860* 1.897$ 

2 

SGDAT DJ=10Q0. * RU=5QQ0. * IG=1®^HT=.5*NT=0$ 

3 

SSORCE I S0=3$ 

7 


1=6 FOR HEADING 
HEADING 

1=1 TO CALL PDATA 
NAMELIST POAT 


1=2 TO CALL GDATA 
NAMELIST GOAT 
1=3 TO CALL GEO 
NAMELIST SORCE 
1=7 TO TERMINATE RUN 
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SAMPLE PROBLEM CARD INPUT FDR SUPER B 


CARO REFERENCE 


S fNPU EM IN= * 09* O‘J= ; XO00a 1> RU :s 5030s»IG = l»IRAN =! «5i?El = 1500®9E2 =2000 » 9 WHT = * 5* 
ETH=g.D0?5«PAAS = «D710^^ALTH=10«89a« NRECPT-50 § 


SAMPLE PROBLEM CARD INPUT FDR STAT 


4 2 

3 l 

L 

1 107021 

1001 A 

5000 * 

2 

5000* 

1000* 


a 05 

0 * 

20* 

0* 

THEP 2 

0 * 

45 * 


F 140 

0 ® 

5* 


2THEP E 

1THFP 

IF 

1 500 e 

2000® 




NAMELIST INPU 


CARD 

REFERENCE 

CARD 

SET 

A 

CARD 

SET 

B 

CARD 

SET 

C 

CARD 

SET 

D 

CARD 

SET 

E 

CARD 

SET 

F 

CARD 

SET 

G 


47 



07/0?/ 7 Q LRC CM SCOPE 1«Q 

15.26*19 *DAPO 47 4» 

15« 26*. 19* 

15*26*19* JOa* I* 100*45000. 
l 5« 26® L 9o J „ LAMB IQTTE 
1 5* 26* 20* RUN ( S) 


6S00:-i3l& 06/02/70 
0702 1515 

LRC COMPUTER COMPLEX 
R0000f90105«2* 
RM 2054 


15*26®?7* LGO* 

15* 26. *41. MEMORY 012600 CM 
15« 26* 43 ® EX l T 

15® 26»44»SP PRINT (OUTPUT ®5) 
l 5® 26® 48<*CPU 0*919780 SEC* 

15.26*48® PPU 29 • 249 536 SEC* 

15® 26* 48. DATE 07/02/70 

17*06*16® DAPQ474. 409 LINES PRINTED® LP21 
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TWOGEV SAMPLE PROBLEM OUTPUT 
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E0= 2000*0 RO = loOC NUMBER OF STEPS- 50 Z= 7 CF=-5„D0 E CJT 0FF=400 o 0 AKK = l c 100 AVG= 4* 18264E+22 


HHfMfMWfMfNJOOlMfri 

OOOOOOOOOO 
++++++++++ 
UJ m uj in ill n i nt HI III in 

oooooooooo 

©OOOOOOOOO 

OOOOOOOOOO 

COvOHCOOsJ-CON'OO 

o*o« 0 eo* 9 « 


f-ir-lC\jr\J(\jr\J(\JC\!C\l<\J 

OOOOOOOOOO 
+ + + + + 4 - + + + + 
UJ LU LU UJ UJ LL! UJ UJ LU UJ 
OOOOOOOOOO 
OOOOOOOOOO 
OOOOOOOOOO 
'tNOhOONvOO^CC 


f-J O OOOOOOc 
OOOOOOOOOO 
+ + 4 * + + +*f- + + + 
UJlLHilLUUjLUlUaiUjUJ 

OOOOOOOOOO 

HHh-(MrnNHOCf>H 

O'h'O'inHvflO'fnNvO 

4 -HO'NHinCMHOO' 

»«•**««••• 


fHrHOOOOOOOO 
OOOOOOOOOO 
++++++++++ 
LU UJ LU LU LIJ UJ LU LU UJ UJ 
OOOOOOOOOO 
o o {Mr-r^rocvcT'rHON 

Orvj'tvj-oj^rnHOO' 

o *n <t- ro <m cm cm r\r ^ 


i— I r— ♦ t — | f\J f\J (\J (\1 Osl CNJ C\J 

OOOOOOOOOO 
+ + + + + + + + + + 
LU UU UJ UJ UJ LL) UJ UJ LLJ UJ 
OOOOOOOOOO 
OOOOOOOOOO 
OOOOOOOOOO 

0°00'OvOO<t‘CO(\]>X) 

0 ©«*Oe*oa» 

HfnO'HM«tlAOcoO' 


rH^HOOOOOOOO 

OOOOOOOOo o 

+ ++ + + + + + + 4 - 
LUUJUJUJUJLUUJUJUJUJ 
OOOOOOOOOO 

CT'Oiror^-comomr- 

* 0 * 0 * 099 * « 



nHHtVcOMNMf'JfMfn f-Hn-nOOOOOOO OO 

CC OOOOOOOOOO o ooooooooooo 

UJ + 4 - + + + + + + + + + + + ++ + + + + + + + 

3 UJUJUJUJUJU-'UJUJUJUJUJ UJLDUJUJUJUJUJUJlJJLUUJ 

o ooooooooooo OOOOOOOOOOO 

Q. OOOOOOOOOOO hNCOO'tflHOOh-M? 

OOOOOOOOOOO (MC'OO'J-UICMJ' O COH 

o o<fO^>tooM'CO'to r-'^'Occunr^-ro'rMoCf'oo 

Z ®O®O 9 O®«® 0 « 

HH 00 m 00 H (M m r. O 00 O' CVJ m C\J (M (Nj M H H 

CL 
Q- 

O O ^ H C.I M M f\' (\l (\J fM (Y' *“ 4 t-*© 0 © 0 © 0 o 00 

1— OOOOOOOOOOO OOOOOOOOOOO 

CD + + + + + + + + + + + + + + + + + + + + + + 

UJUJLUUJLULUUJUJLUUJUJ UJLLJUJUJUJUJUJUJUJUJUJ 

U- CC OOOOOOOOOOO OOOOOOOOOOO 

o uj ooooooooooo wc.'OO'tro'OcoffiHsOo 

3 ooooooooooo (\j in (\j in ro Cvj <\! r-n co O' m 

to O OOO'J W<OC>f C0W>0 C'O^OhtP'tNOO'CO 

Z O Cl e« 99 **«e «*9 <**«©» 

CJ O v£} M r- H W fT| If I vO (J' H ^HOOlO(^MtM(\J(\T-tH 

►H © o 

b- <t Z 

«c <s »-* OHHNC 4 NNM t-*t-<© 00©©0000 

~i <t CL OOOOOOOOOOO OOOOOOOOOOO 

»• 1 + + + + + + 4 + + + + +44 + + t + t + + + 

o o UJLLJUJUJUJUJUJUJUJUJUJ UJ UJ UJ uj IL 1 UJ 111 UJ III O! UJ 

_J K OOOOOOOOOOO OOOOOOOOOOO 

<t co ooooooooooo vo < — > o 'O o > — 1 r~* t~~ r — so • — i 

Oil OOOOOOOOOO^ ^ino*^(^*-Jf^-f r ' 00'0 

<i cc o^up',o^coNoo 4 <t- O) in m O' 04 - w o o oo 

CC CJ a*®®®®®®**® 

O LL 4 CM'DHNn^'Oh U'r-H CC CO.-lU'infOCvlOvJCMCVJfVl,-* 

U_ o UJ 

O UJ 3 

20 — 1 O rH H CVJ CVJ (\J 04 CVJ CM CVJ CO CJ CM-H .-100000000 

U CM <L OOOOOOOOOOO CL ooooooooooo 

#*-*<» o •S- 4 - 4 - 4 - + 4 - 4 > 4 ' 4 - 4 ' 4 - 4 * 4 ' 4 » 4 - 4 * 4 - + + + 4 ' + 

Hr- t/> ujujujujujujujujujujuj cd ujluujlljujujujujujujuj 

* « OOOOOOOOOOO Z O OOOOOOOOOO 

2 - > ooooooooooo r- co r- © *i> o m co f». 

*CL CD OOOroOOOOOO© CL H H (T> N O O H vf) h r 1 O' 

CJ OrT 0<\JO<\ltr>Cvj'00'tCC<M CL 'tHO-OHO“'(MHOCO 

4 t-U.ll UJ *s«eo«» 9 *ee O ® 0 ® a®e«o® 0 e 

2 f^J z (MOviinHHmst-^rccH k Hror-«in* 4 -rn(\jN 04 MH 

* uj to 
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o 


OOOOCJOOOOOOOOOOOOOOOOOOOOOOOOO^— 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

+ + + + + ■+ + + + + + + + + + + + + + + + + + + + + + + + + I I I I t I I I I I I I I I I I I I I 

UJlUUJUjLLJUjLUUJUJUJUJUJUJailtJ|JjUJmUJIJJUJUJUJUJUJLUUJUJUJIJ.JUjUJlijUJUJUJLUCXIUJUjmtJJUJUJUJUJUJUJUJ 
Minj)>fH^ffl(DcoinHincoooo'octnHO^OH(fivt-(Ni^v£)4-o>4-0'OooNOccOHflOin'4'oovON'-iMnooO' 
aCMir*oor-tff'irirc'^iro>tir\I^oo0^oo^^HOOC'cor“iorOr-!OvOs£)rnrfitnmor^'J->4-ruomr-4r-irnr^cjO'vO 
NO^(\iHCJ'r‘um(\jC'CD'0'J-(NJO«)Mrir0H0'v0'd'(Mofl0^rfiHr0O'0r\ihN'OO(fl'0000'oocr-Mf\Hh 
ir> ir, in m u! < ^'j-'4‘4-<4 , mMror<i{^c\jMC\iNNHHHHHOoOoO'^'i-(VOh 'j-wo'>orooooiAHOOint\ioo 
C~< <t9#»#*»»»*>o«»»o®#©0«oeo»o9*eo®#«o®peo®oe©S>o®o«o*o® 

£V OH’^Hr-lr-IHHrlHHHHi-IHHHi-IHHHf-lrlHH^HHHHHCr'a'00'0 0 WCOODhNl^hv£l'0^ininU\4' 

a 
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E0= 2000*0 R0= IsOO NUMBER OF STEPS= 50 1 - 0 CF= 0*00 E CUT 0FF= 2*0 AKK = 0«000 AV0- 



EK 

OS 

PRO 

2® CGOOE+C3 

0. 

0* 

l®9592E+03 

1 «9168E+01 

0. 

1* 9184E + 03 

1<* 9217E+01 

0* 

1.8777E+03 

i « 9268E+0 1 

0* 

U8369E+C3 

1 * 9321E+01 

0® 

1 »7961 E+03 

1 « 937 5 E+01 

0 * 

1* 7553E+ C3 

1® 9431 E+Cl 

0* 

l«7146E+03 

l®9488E+0 1 

0. 

la 67 38E+C3 

1® 9547E+01 

0* 

U6330E+03 

1 ®9607E+Q1 

0® 

la 5922E+C3 

K9669E+01 

0® 

1® 55 15E+03 

1 « 9^325+01 

0. 

1<® 5107E+C3 

1 ® 97 96E + CI 

0® 

1 « 46 99E+03 

1.9863E+01 

0® 

1.4291E+03 

U9931E+01 

0* 

1* 3884E+03 

2® 00C2E + 01 

0® 

1® 3476E+C3 

2 ©0075E+01 

0® 

1 ®3C 68E+03 

2® 01 50E+C1 

Oo 

1 * 26 60E+C 3 

2.0227E+01 

0. 

l«2253E+03 

2®03C6E+01 

0® 

1® 1845E+G3 

2 «0338 E+01 

0. 

1j» 1437E+C3 

2o C478E+C1 

0* 

1« 1029E+C3 

2 ®0568£+0 l 

0® 

l*C622E+03 

2® 0657E+C1 

Oo 

la 0214E+03 

2 ®07 46 E+01 

0. 

9® 80 61E+02 

2o 03 33E + 01 

0* 

9« 39 84E+02 

2 ®0929E+0 1 

c* 

8*9<506E + 02 

2® 1031E+01 

Oo 

8®5829E+02 

2 ®?.133 E+01 

0 » 

8® 17 51E+02 

2®1240E+01 

Oo 

7*7fc73E+02 

2®1345 E+01 

Oa 

7 o 35 96E+02 

2®146CE+C1 

Oo 

6® 95 1 8 E+02 

2 «1 576 E+01 

3® 

6® 5441E+02 

2o 1693E+01 

Oo 

6*1363E+02 

2 ®181 1E+01 

0® 

5*72 86E+02 

2® 1926E+C1 

Oo 

5®32Q8E+02 

2® 2039E+01 

0® 

4®9I 31E+02 

2® 2 1 56E + C 1 

Oo 

4a 5053 E+02 

2®2255E+01 

0* 

4® 09 76E+02 

2® 2348E+Q1 

0® 

3o 68 98E+02 

2 ®2417E+01 

0, 

3® 2820E + 02 

2® 2453E+01 

Oo 

2aP743E*02 

2 ®2437E+0 1 

0® 

2 a 4665E +02 

2® 2329E+01 

Oo 

2® C5 88E+ C2 

2 &2073E+01 

0* 

1® 6510E+02 

2® 1579E+C1 

Oo 

1« 2433E+02 

2 ®0 642 E+0 1 

Oo 

P« 35 51 E+01 

1® 8899E+ 0 1 

0, 

4 a 2776E+01 

1 ®5483E+0 1 

Oo 

2® CCOOE+OO 

8® 0401E + 00 

Oo 
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GDATA 

LENGTH= lOOOaOOO 

MAX NUMBER OF STARTED P ART ICLE S = 500 

HALF SQUARE SIDE= 5000*00 

ATMOSPHERIC CALCULATION 5C0 PARTICLES ISOTROPIC 

BEGIN M *C« AT 06/25/69 17*36*54* 

STARTING RANDOM NUMBERS 
GAUSSIAN *3000 

UNIFORM © 3000 

SOURCE ENERGY INTERVAL IS 1500*0 2000*0 

PERCENTAGE OF NEUTRONS IN SOURCE SPECTRUM IS 40©00 
INITIAL NEUTR* 199 PROT « 301 

STARTING WITH 500 PARTICLES ON RANGE *4GEV TO 2GEV 


NGEN COLL 

P 

N 

PI + 

PIO 

PI- 

1 

499 

301 

199 

0 

0 

0 

2 

542 

1670 

1728 

235 

240 

229 

3 

407 

1543 

1603 

114 

136 

124 

4 

243 

1C 52 

1102 

7C 

69 

75 

5 

1C6 

5E4 

651 

21 

30 

42 

6 

42 

259 

259 

8 

12 

14 

7 

10 

ICO 

99 

1 

1 

3 

8 

3 

27 

29 

1 

0 

l 

9 

0 

7 

11 

C 

0 

1 

END 

M* C « AT 

06/25/69 17 

• 42,05, 



NO® 

OF LOW 

NEUTRONS 

AND 

PROTONS 

IS 

9063 

NO* 

OF PROTONS SLOWED TO 400*0 

MEV 

IS 255 


END OF RUN, NO* OF RECORDS IS 


25219 



06/25/69 LRC CM SCOPE 9 0 C 6600C-121K C4/16/69 
17, 31* 16. CAPS 7, 56. - 0625 1504 

17o31ol6« L RC COMPUTER COMPLEX 

17. 31 cl 6, JOB. 1.400, 72000. ROOOO , 3 1 67 1 , 7, 

37. 2UI6.JAY J. LAMBTOTTP, RM 2054 

17,31.36. RUN ( S } 

17. 31« 4 7eP0LLQUT WAITING ON ALL QUIET 
17. 31.48. ROLLOUT INITIATED. 

17. 31 .50. ROLLOUT COMPLETED. 
i7.3?,48«RCLL IN INITIATE!). 

17. 33.02. POLLIN COMPLETED. 

17. 33.02. REQUEST 0UM.M2.HY* 212054. ROL.JJI , MC 2GEV 

17.23.02. AIR 
17,33.92. (65 ASSIGNED! 

17. 33,02. COPYRE (DUMM2 » TAPE C ?1) 

17, 33. C6. REWIND (TAPE9) 

17 0 33o07«DR0PFIL< DUMM 2 ) 

17. 33.10. RELEASED DUMM2 

17. 33. 14. ROLLOUT WAITING ON ALL QUIET 

17,33, 14, ROLLOUT INITIATED* 

17® 33o 1 6. ROLLOUT CCMPLE TE D. 

17. 34.56. RCLLIN INITIATED. 

17* 36.11, POLLIN COMPLETED. 

17. 36. 12, REQUEST T A p 24 , HY * 5090 1 1 , RI L, J JL , TWOG F V 
17. 26. 12. TEST CASE2 

17.36.12, (61 ASSIGNED) 

17. 36o 1 2 oL GO. 

17. 36.46.MEM0RY 065100 CM 

17® *2. 05. STOP 

17. 42.0 5 *SPPR I NT (OUT PUT, 5) 

17® 42. 07. UNLOAD (TAPE4) 

17. 42. 13, CPU 112.691689 SEC » 

17. 42* 1 3. PP U 224. 813056 SEC * 

17* 42® 13® CATE 06/25/69 

18® 39® 3 5® DAP5456® PRINT-PP C2032 LINES, LP 25 
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PROTOS SAMPLE PROBLEM OUTPUT 



NEW PHYSICAL DATA 


* CM CM CM CM < 

) O O O O I 

+ + + + 

I LU UJ UJ UJ I 

> o o o o < 

f o o o o < 

> © O O O < 

IHCOO-tl 


ir-HOOOOOOOOO 
•OOOOOOOOOO 
+ + + + + + + + + + 
lUIUJUJUJUJUJUJUfUJUJ 
IOOOOOOOOOO 
irHf^(Mror*-^-«oo s i-Hro 
r-O'^H'OafntM'OfO 
HO'(MHlfU\JHOO'<D 


H H N CM ty (\| CNJ N W CO (f) 
OOOOOOOOOOO 
+ + 4 , + + , f + + + + + 

UJUJLUUJUJIJJUJIUUJUJUJ 
OOOOOOOOOOO 
OOOOOOOOOOO 
OOOOOOOOOOO 
NOh® (Vl'OO'tffla) 
• ••••««•««« 


>00000 
) O O o o o 
' 4- + + 4 + 

J UJ UJ UJ UJ UJ 

>00000 

JMO'HCf N 
4 CM > 4 " <r CM 
) m mo cr a 


MHMNWW 

000000 

4- + + 4- 4- 4- 
UJ UJ UJ UJ UJ UJ 
OOOOOO 

000000 

OOOOOO 

OCOO'O'OO 


CM CM CM CM rO 

OOOOO 
4- 4 + 4* 4- 
UJ UJ UJ UJ UJ 
OOOOO 
OOOOO 
OOOOO 
4 CON ^ 4 


MfOO'HlM« 4 -l^»OOOC'CM 


IOOOOOOOOOO 
'COOOOOOOOO 
+ + + + + + + + + + 
IUJUJUJUJUIUIUJUIUJUJ 
IOOOOOOOOOO 
is 4 - 0 'ooorO<fC 7 'cacovO 
vOioor-ooiTivOini^o 
ff) O' vO (TW) m H o O' CD 
>•••«• ■•«••« 
I M MO 4 rONfM(\INHH 


Or-IOCMCMCMCMCM(MCMrO 
OOOOOOOOOOO 
+ + + + + 4‘ + + + + 4- 

UIUIUIUIUIUIUJUJIJUUJUJ 

OOOOOOOOOOO 
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oa'oor^-4‘r-tr^cvjr“Omif\f^oir'roovo^inr*oooor->4 , o>3rr-<jc?ootno'Topoo'0^iintom<\ja'a'4-c«o<o 
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a crnf\J(vJ^(NJogcvfM^rg(^(Nirgrsjf\JCsjfvifycjf\if> 4 Cv!cviCMc^c\i^ 4 r- 4 r 4 T~if- 4 ^iHrH^ < -ip-if^^Hwcr'oor**»oinNj' 
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1 « 7286 E-HD 1 3 , 7497 E - 01 8 .A 776 E -03 

9 * 1 A 29 E 4 -Q 0 2 . 5435 E -01 A. 8023 E -03 

1 « 00 C 0 E-S -00 1 . 1661 E -01 1 . 7279 E -03 



GCATA 


LENGTH 1000*000 

MAX NUMBER OF STARTED PARTICLES^ 0 

HALF SQUARE SIDE= 5000.00 
3 

ATMOSPHERIC CALCULATION TEST CASE 
BEGIN M*C.» AT 06/26/69 19*38*42* 

STARTING RANDOM NUMBERS 
GAUSSIAN *3000 

UNIFORM * 3000 


SOURCE PARTICLES FROM BIOGRAPHY TAPE 
SPLIT 2-1 ON FIRST NUCLEAR INTERACTION 


GEN 

COLL 

P 

N 

1 

0 

0 

0 

2 

1490 

1394 

1430 

3 

2591 

3096 

3675 

4 

2965 

4CC1 

5037 

c 

2751 

3818 

5205 

6 

2082 

3112 

4543 

7 

1330 

2229 

3240 

8 

718 

1258 

2C23 

9 

366 

644 

1C78 

10 

167 

313 

53C 

11 

66 

145 

242 

12 

26 

47 

96 

13 

li 

21 

37 

14 

3 

6 

15 

15 

2 

1 

3 

16 

0 

1 

2 


END M*C* AT 06/26/69 19.5U49* 

NUMBER Of LOU NEUTRONS IS 15809 

END OF RUN, NO* OF RECORDS IS 182333 
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06/26/69 LRC CM SCOPE 3,0 6600B— 65K 04/16/69 

19,C8*C9,CAP5577* - 0626 1455 

19.C8.09. LRC COMPUTER COMPLEX 

19*C 8,C 9, JCB* 1*2000,63000* ROOCO *31671 « 2 , 

19«C8eC9*JAY J, LAMBIOTTE* RM 2054 

19*08, 09«RUN(S I 
19®C8,27*LINECNTI 10000) 

19 *C 8. 29, ROLLOUT WAITING ON ALL QUIET 
19*C8«30. ROLLOUT INITIATED, 

19*08,31. ROLLOUT COMPLETED. 

19,14, 55.RCLLIN INITIATED* 

19*14, 56.R0LLIN COMPLETED* 

19*1A,56*RE0UEST DUMM?,HY«* 1202019, ROL* JJL, MC DAT 

19 , 14 * 56« A AIR 

19*1 A® 56* (67 ASSIGNED 

19*14* 57*C0PYBFI0llMM2«TAPE9«lJ 

19,14,59*REWIN0 ( TAP E9) 

19*15*00«CRDPE IL (0UMM2 ) 

19, 15*02* RELEASED DUMM2 
19, 15®C2,SET INDE , 

19*15, 04 8 R CL LOUT WAITING ON ALL QUIET 
19*15,05, ROLLOUT INITIATED* 

19, 1 5, 07® ROLLOUT COMPLETED * 

19,18, 43 <»R CL LIN INITIATED, 

19*18, 45, ROLLIN COMPLETED, 

19.16.45, REQUEST TAPE7?HY, 509011 ,ROL*JJL,TWOGEV 
19,18 ,45®TES T CASE 

19.18.45, 167 ASSIGNEOS 

19, 16, 4 9* ROLL OUT WAITING ON ALL QUIET 
19, ie, 50, ROLLOUT INITIATED, 

19.18.51, ROLLOUT COMPLETED, 

19.22.51, ** 211050 **, 

19.26. C6.RCLLIN INITIATED, 

19.26, C8»P CL LTN COMPLETED, 

15,26,11® 

15,38,2l»REQUEST TAPE4* HY, SAVTP ,RI L * J JL , PROT TE 
15,36* 22® ST CASE2 
15*38,22, (64 ASSIGNED) 

19.38.22, LGO , 

19®38®39,M EMORY 056000 CM 
19, 39, 03, MT 64 WPE RECOVERED, 

19.39. 23, MT 64 WPE RECOVERED, 

19,39, 47.MT 64 WPE RECOVERED® 

15*41. 10. MT 64 WPE RECOVERED® 

19,43, 10®MT 64 WPE RECOVERED® 

19 *4 3, 27* MT 64 WPE RECOVERED, 

19, 51*49, STOP 

19,51 ,50® SPPRI NT! OUTPUT *5) 

19.51*51, UNL0ADITAPE4J 
19«52*26*UNL0AD(TAPE7) 

19,52, 27, CPU 591,370217 SEC® 

19*52*27, PPU 1083*629568 SEC, 

19, 52* 27, CATE 06/26/69 

2C.11.01* CAP5577. PRINT-PP 01932 LINES, LP 21 
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SUPER B SAMPLE PROBLEM OUTPUT 
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LOW NEUTRONS 


EMIN ,090 DU 1000*000 RU 5000*000 IG 1 

ETH 2®5000E-03 PAAS 7*1040E~02 ALTH U0898E+01 

INC® ENERGY INTERVAL IS 1500,000 2000,000 WHT IS ©50 


RANDOM START = ®50000 


X 

Y 

Z 

THETA 

PHI 

E 

IK 

SOURCE 

WT 

GEN 

0 ® 

0 ® 

0 ® 

2» 386 1E+01 

-U8401E+01 

1® 6000E+03 

4 

1 

1 o 000 

1 

5* 5432 E+00 

-1® 8668E+00 

1 ® 2960E+0 1 

2 © 421 3E+01 

-1 « 8710E+01 

1 • 5755E+03 

1 

1 

1 a 000 

1 

1« 2896E+01 

-4* 4084E + 00 

3© 0341E+01 

2 ® 418 8E+01 

-1 ® 9270E+01 

1© 5429E+0 3 

1 

1 

uooo 

1 

1*6778E+01 

— 5© 7288E+00 

3 ® 9400E+01 

2 ©4311E+Q1 

-1 © 8144E+01 

1 * 5254E+C3 

1 

1 

1 © 000 

1 

1.6778E+01 

-5*7288E+00 

3 • 9400E+01 

2 ©431 1E+01 

-1 © 8144E+01 

1 « 5254E+03 

3 

1 

1 © 000 

1 

1© 6778E+01 

—5® 7288E+00 

3 « 9400E+0 1 

1 ® 2354E+02 

7»6609E+01 

6»9877E+0 1 

-14 

1 

1®000 

2 

1*8765 E+01 

2® 6152E+00 

3 ©3713E+0 1 

1 ® 2354E+02 

7® 66C9E+01 

4® 2775E+01 

-11 

1 

1 a 000 

2 

2© 0317E+01 

9® 1340E+00 

2® 9270E+01 

1 © 2354E+02 

7®6609E+01 

2 ©OOOOE+CO 

-11 

1 

1 a 000 

2 

1 o 6778 E+01 

-5© 7288E+00 

3 «9400E+01 

3® 1877E+01 

7* 2714E+01 

1© 530 5E+02 

0 

1 

1 © 000 

2 

1® 6778 E+01 

— 5® 7288E+00 

3 • 9400E+01 

2* 1429E+01 

-5 © 7294E+01 

1 ® 3207E+02 

14 

1 

la 000 

2 

la 7552E+01 

— 6 ® 9337E+00 

4®3048E+01 

2© 1429E+01 

-5®7294E+01 

1 a 2433 E+02 

11 

1 

1 a 000 

2 

2® 1283E+01 

-U2744E + C1 

6 © 0641 E+01 

2© 1429E+01 

-5«7294E+01 

8 ® 3551E+01 

11 

1 

1 a 000 

2 

2«4340E+01 

-1® 7504E+01 

7 ® 50 54E+01 

2® 1429E+01 

-5«7294E+01 

4* 2775E+01 

11 

1 

1 « 000 

2 

2® 5927E+01 

“1 © 9976E+01 

8 ®25 38E+0 1 

2© 1429E+01 

-5®7294E+01 

2 © OOOOE + OO 

11 

1 

1 e 000 

2 

1® 6778 E+01 

-5# 7288E+00 

3® 9400 E+01 

9® 1867E+01 

6 ® 81 47 E+01 

2® 178 5E+01 

-4 

1 

1 a 000 

2 

4, 6104 E+01 

6 ® 7396E + 01 

3 ®6831 E+01 

9© 1867E+01 

6 ® 8 147E+01 

2© 1785E+01 

-1 

1 

1 « 000 

2 

4«6104E+01 

6 ® 7396E + 01 

3 ©6831E+01 

9 e 1867E+01 

6 ® 8 147 E+01 

2 « 1785E+01 

-3 

1 

1 a 000 

2 

4® 6104E+01 

6 « 7396E + 01 

3 ©68 31E+01 

9® 7754E + 01 

9» 9798E+01 

1 © 8427E+01 

-4 

1 

a 500 

3 

1® 9520E+01 

2® 2133E + 02 

1 ® 55 60E+0 1 

9 ®7754E + 0 1 

9® 9798E+0I 

1 © 8427E+01 

-1 

1 

® 500 

3 

1® 9520E+01 

2® 2133E+02 

1 ® 55 60E+01 

9® 7754E+01 

9* 9798E+01 

1 © 8427E+01 

• - 3 

1 

a 500 

3 

lo 9520E+01 

2® 2133E + 02 

1 »55 60E+0 1 

1 ® 0843E+02 

-1 © 5893E+02 

1 © 176CE+00 

4 

1 

*500 

4 

1«9518E+01 

2® 21 33E+02 

1 * 5559E+0 1 

1 * 1354F+G2 

-1®6371 E+02 

laOOOOE+OO 

1 

1 

a 500 

4 

4® 61 04E+0 i 

6®7396E+01 

3 ©6831E+01 

7 © 6090E+01 

6 ®9852E+01 

1 ®6410E+01 

-4 

1 

1 a 000 

3 

4«771 6E+01 

7® 1788E+01 

3® 7990 E+01 

7 « 6090E+01 

6 ®9852E+01 

1* 6410E+0 1 

-1 

1 

1 © 000 

3 

4®7716E+01 

7® 1788E+01 

3 «7990E+01 

7® 6090E+01 

6«9852E+01 

1 ® 6410E+01 

-3 

1 

1*000 

3 

4©7716E+01 

7® 1788E+01 

3 «7990E+01 

8© 174 2E+01 

6 ® 9793E+01 

1 © 6399E+01 

-4 

1 

1 a 000 

3 

5*0855 E+01 

8 ® 0316E+01 

3 ®93C9E+01 

8© 1742E+01 

6 © 9793E+01 

1 © 6395E+01 

-1 

1 

1 *000 

3 

5* 0855E+01 

8 © 0316E + 01 

3 ®9309E+01 

8 ® 1742E+01 

6 ® 9793E+01 

1 © 6399E+01 

-3 

1 

1 a 000 

3 

5® 0855E+01 

8 ® 0316E+01 

3 ©9309E+01 

1 « C719E+C2 

6 © 5582E+01 

4©4907E+00 

-4 

1 

1 a 000 

3 

5® 8663E+01 

9® 7514E+01 

3 ® 3464E+0 1 

1 a 0719E + 02 

6 # 5582E+01 

4*4907E+00 

-1 

1 

1 e 000 

3 

5# 8663 E+01 

9® 7514E + 01 

3« 3464E+C1 

1 © 0719E+02 

6 ® 5582E+01 

4®4907E+00 

-3 

1 

1 a 0 00 

3 

5* 8663 E+01 

9© 7514E+01 

3®3464E+01 

1 a 4525E+02 

2« 2456E+01 

4« 2657E+00 

-4 

1 

1 a 000 

3 

6 ® 7675E+01 

1® 01248+02 

1 ® 9408E+0 1 

1 ©4525E+02 

2® 2456E+01 

4© 2657E+00 

-1 

1 

la 000 

3 

6 * 7 675 E+01 

1® 0124E+02 

1 ©94C8E+01 

1 ©4525E+02 

2« 2456E+01 

4® 2657E+00 

-3 

1 

1®000 

3 

6 ® 7675E+01 

1 « 0124E+02 

1 ©9408E+0 1 

1 ©6331E+02 

-1 ® 4702E+02 

4® 0708E+00 

-4 

1 

1 a 000 

3 

5© 9101 E+01 

9® 5676E + 01 

- 1 ©4679E+01 

1 « 6331E+02 

-1 ©4702E+02 

4® 0708E+00 

-2 

1 

1 a 000 

3 

1* 6778 E+01 

-5®7288E+00 

3 ©9400 E+01 

3©9889E+01 

-3© 5162E+01 

1© 1778E+03 

-4 

1 

1*000 

2 

4© 052 8 E+01 

— 2® 2459E+01 

7®41 57E+01 

3* 9889E+01 

-3© 5162E+01 

1© 1778E+03 

-1 

1 

1 a 000 

2 

4® 0528 E+01 

-2® 2459E + 01 

7«41 57E+01 

3«9885E+01 

— 3©5162E+01 

1 ® 1778E+03 

-3 

1 

1 a 00 0 

2 

4®C528E+01 

“2® 2459E+01 

7*41 57E+0 1 

8 * 8604E+01 

2 «3897E+01 

4®7486E+01 

14 

1 

1 ® 000 

3 

4® 2163E+01 

-2® 1734E+01 

7©42 C1E+01 

8 ® 8604E + 01 

2 ® 3897E+01 

4® 2775E+01 

11 

1 

1 ® 0 00 

3 

4®9512E+01 

—1® 8478E + 01 

7®4396E+01 

8 « 8604E+01 

2 ® 3897E+01 

2 © OOOOE+OO 

11 

1 

1 a 000 

3 

4® 0528 E+01 

-2© 2459E+01 

7 ®4157E+0 1 

4* 5658E+01 

-4® 2576E+01 

4* 7080E+02 

-4 

1 

1 © 00 0 

3 

5® 0397E+0 1 

-3® 1526E+01 

8®7254E+01 

4© 5658E+01 

-4© 2576E+01 

4®7080E+02 

-1 

1 

1 a 000 

3 

5<* 0397E+01 

-3® 1526E+01 

8 ©7254E+01 

4® 565 8 E + 01 

— 4« 2576E+01 

4® 7C80E+02 

-3 

1 

1*000 

3 

5® 0397E+01 

-3® 1526E+01 

8*7254 E+01 

6 © 262 1E+01 

2 ® 5741E+01 

1 ® 6882E+02 

-4 

1 

1*000 

4 

2©0760E+02 

4® 4270E+01 

U7763E+02 

6 ® 262 1E+G1 

2® 5741E+01 

1 ® 6882E+02 

-1 

1 

1*000 

4 

2® 0760E+02 

4© 4270E + 01 

1 ©7763E+0 2 

6 © 262 1E+C1 

2© 5741E+01 

1 ©6882E+02 

-3 

1 

1®000 

4 

2*0760E+02 

4® 4270E+01 

1 *77 63 E+02 

1 * 2997E+02 

-1 *2551 E+02 

8®7780E+00 

-4 

1 

1©000 

? 

2*C541E+02 

4* 1205E+01 

1 * 7448E+02 

1 * 2997E+02 

—1® 2551E+02 

8 ® 778CE+00 

-1 

1 

1«000 

c 


END OF RUN? NO© OF RECORDS IS 400701 
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06/27/69 LRC CM SCO PE 3*0 6600D-131K C4/16/69 

10* 21® 03* DAP 1603® 

10® 21® 03® LRC COMPUTER CCMPLEX 

10* 21® 03* JOB a * 600 * 47000* ROOOO * 31671 , 2 , 

10® 21*03® J® LAMBIOTTE RM 2054 

10* 21® 04* RUN { S ) 

10®21®18<*SETINDF» 

10®21«19® L INECNT ( 10000 ) 

10*22© 09® REQUEST TAPPli *HY © 114018 *ROL * JJL ? NUC DA 
10® 22® 09*TA (400 MEV DOWN I 
10® 22® 09® (55 ASSIGNED) 

10© 22® 10*C0P YBF ( TAPP11? DI S ) 

10® 22«44*C0PYBR( TAPP 11 *DI S »24 ) 

10®23®23®C0PYBR(TAPPllfTAPE10#2) 

10® 25* 01« TAPE 107021 ® ® »© ® 

10® 25* 21® REQUEST TAP E3? HY © SAVTP 9 R!LtJJL ? SB TEST 
10*25*21® CASE2 
10*25*21® (65 ASSIGNED) 

10*26*07«REQUEST TAPE4*HY® 21 1050? ROL? J JL? PROT TE 
10*26® 07* ST CASE 
10*26® 07® (66 ASSIGNED) 

10® 26* 07® LGG® 

10® 26® 21® MEMORY 036500 CM 
10® 26® 27®MT 65 WPE RECOVERED® 

10*37o47*MT 65 WPE RECOVERED® 

10*37®48®MT 65 WPE RECOVERED® 

10®37® 50®MT 65 WPE RECOVERED® 

10® 37* 51® MT 65 WPE RECOVERED® 

10® 37® 52® MT 65 WPE RECOVERED® 

10* 37* 53* MT 65 EOT 

10*42*04® TAPE 107021 CONT TO 111018® « ® ® ®« ® «» ® © ® *® * 

10®42®33« 

10® 44* 34® EXIT 

10® 44® 34® EOF ENCOUNTERED BY NAMELIST 
10* 44® 34* SPPR I NT (OUT PUT 9 5 ) 

1 0* 44* 3 5® UNLOAD ( TAP E3 ) 

10*44® 36® UNLOAD! TAP E4) 

10* 44* 37® CPU 290® 765124 SEC® 

10*44® 37»PPU 1310*113792 SEC® 

10»44® 37©DATE 06/27/69 

11*04*25® DAP1603* PRINT-PP 01300 LINES,LP 26 
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STAT SAMPLE PROBLEM OUTPUT 


63 



4 PLANES 


2 VARIABLES 


3 STATISTICS 


OPTION 1 


WEIGHTING FACTOR 2 


STATISTICS REQUESTED FOR NEUTRONS 
STATISTICS REQUESTED FOR PI + 

4 BINS/WORD AND 14 BITS/BIN 
THE MAXIMUM VALUE FOR ONE BIN IS 


VARIABLES USING 
THEP E 

NO® LOCATIONS NEEDED PER PLANE IS 140 
4 PLANES USED IN THIS PASS 


8191 «5 
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THE INCIDENT ENERGY INTERVAL FOR THIS CASE IS 1500*00 2000*00 

IK= 1C000 NO* RECORDS USED= 272993 NO* LOW N EUTRONS= 90613 

THE KO * CF TINES OVERFLOW NOTED IS 0 
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PLANE LOCATED AT Z= 50*00 
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PLANE LOCATED AT Z = 50.00 


STATISTIC 

VARIABLE 

N(THEP* 

) 

FOR NEUTRONS 
START DELTA 

NO. 

BINS 


THEP 


C«OG 

45.0000 


2 

THEP 


717,50 

865.50 





STATISTIC 

VARIABLE 

NITHEP, 

1 

FOR PI* 
START 

DELTA 

NO. 

BINS 


THEP 


0*00 

45® OOOO 


2 

THEP 


31 ® 00 25.00 
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PLANE LOCATED AT Z= 5D o 00 
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PLANE LOCATED AT 2 = 100*00 
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STATISTIC NCTHEP.E 
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PLANE LOCATED AT 1 = 100«00 


STATISTIC 

N (THEP * ) 






FOR NEUTRONS 



VARIABLE 


START 

DELTA 

NO* BINS 


THEP 


o 

o 
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45» Q0C0 
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THEP 
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STATISTIC 

N (THEP ? 1 

► 

FOR PI+ 




VARIABLE 


START 

DELTA 

NO* BINS 


THEP 


0*00 

A 5* 00C0 

2 

THEP 


35*00 

29*00 
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PLANE LOCATED AT Z= 100. OQ 
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06/27/69 LRC CM SCOPE 3«0 66 CCB — 65K 04/16/69 

15©55®28®DAP5679« - 0627 1334 

15© 55.® 2 8© LRC COMPUTER COMPLEX 

15® 55®28®JOB,1,400 ,65000© ROOOO *3X671 * 2 * 

15* 55« 28# J® LAMBIOTTE RM 2054 

15* 55® 29© RUNf SI 

15« 56*06© ROLLOUT WAITING ON ALL QUIET 
15© 56® 08« ROLLOUT INITIATED© 

15® 56® C9®R0LL0UT COMPLETED® 

16*01®04®R0LLIN INITIATED® 

16«C1«06« ROLL! N COMPLETED® 

16« 01 *06 * REQUEST TAPE4, HY© 107021,RC, JJL,TEST CASE 
16® 01® 06® 

16® 01® 0 6* (61 ASSIGNED S 
16© 01® 0 6 « LGC® 

16® C 1® 27® MEMORY C63200 CM 
1 6® 06® 24 ® MT 61 EOT 
16© C9® 34© EXIT 
16® C9*35©$PPRINT( OUTPUT, 5) 
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RESULTS AND DISCUSSION 


The calculations performed by these programs have been compared with other sim- 
ilar calculations and with experimental data. These comparisons differ in two fundamen- 
tal ways. The calculations with which comparisons are made used the same nuclear 
interaction data as PROPER 3C, namely, those generated by Bertini; thus the accuracy of 
the transport calculation itself is checked. The comparison with experimental data is an 
overall check on both the transport calculation and the nuclear interaction data. 

The radiation dose in the atmosphere for a rigidity spectrum with rigidity param- 
eter P 0 of 100 (where c is speed of light) was calculated by Leimdorfer, Alsmiller, 

V 

and Boughner (ref. 9) and is compared in figure 5 with the results obtained from 
PROPER 3C. The absorbed dose (rad) should be used for comparison since the calcula- 
tion of the biological dose (rem) is model dependent. Further comparisons with results 
of Armstrong, Alsmiller, and Barish (ref. 10) for the February 1956 solar proton event 
are shown in figure 6. The upper-limit spectrum for this event was different for the two 
calculations. Additional details of the results of these comparisons are given in 
reference 11. 

The resultant dose rates in the upper atmosphere at high latitudes due to galactic 
cosmic rays measured 1 year after solar minimum of cycle 20 are compared with the 
calculations of PROPER 3C in figure 7. In order to accomplish this calculation, assump- 
tions had to be made for the interaction of the heavy primary particles. Also, the incident 
spectra are not precisely known (especially for energies less than a few hundred MeV and 
greater than a few GeV). A detailed discussion of the dose rates and assumptions is given 
in reference 12. 

Although PROPER 3C has been and can be used to obtain useful results, it does have 
some restrictions which should be restated. Several geometric restrictions of the code 
may be overcome with modifications. First, the geometry subroutine STATI assumes 
either a circular or square cross section for the slab. A more complicated geometry 
subroutine package could be inserted in its place. Second, the code is written for a single 
homogeneous slab of material as a shield. To include calculations for multilayer com- 
posite slabs, it would be necessary to transport all the particles through the first slab 
and then to use the output biography tape as a source tape on the next slab. This proce- 
dure can be accomplished with minor programing changes. 

At the time of this publication, the data available for the high-energy portion of the 
code (program TWOGEV) are not substantial. The data which are available for oxygen 
and aluminum do not include the pion interaction cross sections. It is planned to update 
the data library whenever possible. Even though the code does not depend upon the con- 
tent of the cross-section data, it is highly dependent upon the form in which the data have 
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been tabulated. This dependence complicates the introduction of data from sources other 
than those referred to herein. 

Finally, several features of the code make it machine dependent (CDC series 6000 
computer systems). Nearly all these features are the result of attempts either to reduce 
core storage requirements or to reduce the number of output biography tapes. To reduce 
core storage requirements, the FACT array is packed six values to a word for TWOGEV 
and four values to a word for PROTOS. The packing scheme is highly dependent upon the 
word length (60 bits) of the CDC computer. STAT also contains word packing with the 
same dependence. The data to be written on the biography tape are first packed two val- 
ues to a word by a COMPASS subroutine DPACK and then written on tape by means of a 
blocking subroutine RECOUT both of which are machine dependent. Subroutines RECIN 
and UNPACK perform the inverse operations to read the data. 

CONCLUDING REMARKS 


The PROPER 3C transport code has been successfully used to evaluate the differ- 
ential spectrum of secondary particles emerging from a large number of materials when 
these materials are bombarded by energetic nucleons of an arbitrary source configuration 
The differential spectra after multiplication with appropriate quality factors are inte- 
grated to give the radiation dose deposited at any given point in space. 

The programs in PROPER 3C are flexible enough to permit many types of radiation 
transport calculations with little or no modifications by the users. The programs have 
some limitations, among which are the simplified geometry of the transport medium and 
the lack of a detailed pion transport. However, both limitations can be removed without 
much difficulty. 

Because of the block structure of the program, one can introduce any degree of 
geometric complexity by changing subroutine STATI. 

The mechanism for pion transport is included in the program but the program does 
not simulate the transport of pions because the pion-nucleus cross sections are not 
tabulated. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 27, 1971. 
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APPENDIX A 


PROGRAM ABSTRACTS 


The following program abstracts give a brief description of the development of each 
of the respective programs in PROPER 3C. Copies of these abstracts are kept on file. 
Any communication concerning any of these programs should use the program number 
(e.g. R2090) as a means of identification. 
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APPENDIX A 


01 4 

01 7 PROGRAM NO 

LAR 

R2090 


NASA- LANGLEY RESEARCH CENTER 


COMPUTER PROGRAM ABSTRACT 


101569 


01 20 TITLE OF PROGRAM (61 CHARACTERS MAXIMUM) 

Preparation of Monte Carlo Data Tape for R2Q92 




02 27 

LANGUAGE 
NO. 1 

02 '32 

LANGUAGE 
NO. 2 

FeSR6M 

■ 



High-Energy Hucleon-Pion Transport^ Monte Carlo 


WHO TO CONTACT ABOUT THE PROGRAM 

j 05 28 SITE JOS 31 ORGN CODE |CS 39 PROJECT NO I OS 45 


J. J. Lambiotte LAR H.I50 BIP32k 


06 48 STATUS 

1 1 a. under development 

l I B. operational 
S3 c. completed 


05 49 

a A. THIS PROGRAM 

is not for 

SHARING 



05 58 REVISION CODE 

I 1 A REVISION 
C3 B CANCELLATION 



ABSTRACT 

This progr a m is a modification of R2089 and is use d to prepare! 
a tape for E2092 which is a transport program requiring 1 

nuclear data for pions in addition to neutrons and pro tons , 

These data include energies extended to 2000 MeV * . 
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01 4 

01 7 PROGRAM NO 

LAR 

R2089 


NASA- LANGLEY RESEARCH CENTER 


COMPUTER PROGRAM ABSTRACT 


01 14 DATE 

IOI569 


01 20 TITLE OF PROGRAM (61 CHARACTERS MAXIMUM* 

Preparation of Monte Carlo Data Tape for B2091 


02 32 02 37 KEY WORDS *8 MAXIMUM . SEPARATED BY COMMAS) 

™ E Nucleon Transport, Monte Carlo 


02 14 

CATEGORY 

K 



05-'48 STATUS 

05 49 

1 1 A. UNDER DEVELOPMENT 

CZ3 B. OPERATIONAL 
1X3 C. COMPLETED 

\ 1 A. THIS PROGRAM 
IS NOT FOR 
SHARING 



ABSTRACT 


This program prepares a Monte Carlo data tape which contains 
nuclear data for a particular transport medium. The data are 
for protons a nd neu trons -with energies of 400 MeV (or less). 
This tape is input to R2091* 
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NASA -LANG LEY RESEARCH CENTER 


01 4 01 7 PROGRAM NO 

LAR R2092 


01 14 DATE 

101569 


01 23 TITLE OF PROGRAM (61 CHARACTERS MAXIMUM) 


2000 MeY Monte Carlo Nucleon Transport 


02 26 

CATEGORY 

02 27 

LANGUAGE 
NO. 1 

BO 

F0R6M 


02 37 key WORDS <8 MAXIMUM. SEPARATED BY COMMAS) 

Monte Carlo, Nucleon-Pi on Transport, High-Energy Transport 


05 58 REVISION CODE 
CD A REVISION 
CD B CANCELLATION 


05 -48 STATUS 

05 49 

1 1 A. UNDER DEVELOPMENT 

1 1 B. OPERATIONAL 

Pyp C. COMPLETED 

[ 1 A. THIS PROGRAM 

IS NOT FOR 
SHARING 


ABSTRACT 


Progra m B20 9 1 has been modified to alloy energie s u p to ! 

2000 MeV and to inclu de the transport of pions. The nucleons 

are tjreated only ahove~¥00 MeY„ Any particles created or 

slowed to ielowkOO MeY are treated under R2091* 1 
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01 4 

01 7 PROGRAM NO 

LAR 

R2091 


NASA- LANGLEY RESEARCH CENTER 


COMPUTER PROGRAM ABSTRACT 


101569 


01 20 TITLE OF PROGRAM |6] CHARACTERS MAXIMUM] 


02 27 

LANGUAGE 
NO. 1 


F$R6M 




400 MeY Monte Carlo Nucleon Transport 


Electromagnetic Interactions 


WHO TO CONTACT ABOUT THE PROGRAM 

1 05 28 SITE (05 31 OfiGN CODE }C5 39 PROJECT MX j 05 45 


05 14 CONTACT 05 28 SITE 05 31 ORGN CODE 05 39 PROJECT h 

J. J. LamMotte LAR 11.150 RBP52lj. 



05/48 STATUS 

05 49 

1 1 A. UNOER DEVELOPMENT 

CZ3 B. OPERATIONAL 
f^*~l C. COMPLETED 

| \ A. THIS PROGRAM 
IS NOT FOR 
SHARING 
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01 •! 01 7 PROGRAM NO. 

L,AR R2093 


NASA- LANGLEY RESEARCH CENTER 


COMPUTER PROGRAM ABSTRACT 


01 20 TITLE OF PROGRAM (61 CHARACTERS MAXIMUM) 


Monte Carlo Transport of Low Energy Neutrons 

02 32 02 37 KEY WORDS -8 MAXIMUM . SEPARATED BY COMMAS) 

LANGUAGE 

2 Low-Energy Neutron Transnoi 


101569 


PARENT PROGRAM 


K 1 MSC 


CPS60 


WHO 

TO CONTACT 

ABOUT THE PROGRAM 


05 14 CONTACT 

05 28 SITE 

05 31 ORGN CODE 

C5 39 PROJECT NO. 

05 45 
NASA 
CENTER 

J. J, Lambiotte 

LAR 

11.150 

KEfP^ 


DATES 1 

05 58 REVISION CODE 

05 50 INITIATED 

OI67 

05 54 COMPLETED 

1069 

1 \ A REVISION 

CZJ B CANCELLATION 


05 48 STATUS 

05 49 

1 1 A. UNDER DEVELOPMENT 

1 1 B OPERATIONAL 

£□ C. COMPLETED 

( | A. THIS PROGRAM 
IS NOT FOR 
SHARING 


This program uses 
transport of ener 


(gx_ j.Qlx g las 1 4 g ± - Inela sti c ? an d_f i^s lm_l.nt^ ,.are-_us 

h y the .pro gram » Ba Q93-d.s_.the .final ...link . in. the transport 

initiated in B2092 % it outputs a final history tape which is 
a nalyzed by B209^« __ 

















LAR R2094 


APPSNDXX A 

NASA-LANGLEY RESEARCH CENTER 
COMPUTER PROGRAM ABSTRACT 


02 26 02 27 02 32 02 37 KEY WORDS (8 MAXIMUM . SEPARATED BY COMMAS) 

CATEGORY LANGUAGE LANGUAGE , _ . , , ___ 

no. 1 no. 2 Monte Carla Transport , Particle FI 


Fjz 5R6M CPS60r 


01 14 DATE 

101569 


l - 1 — — — — 

01 20 TITLE OF PROGRAM (61 CHARACTERS MAXIMUM) 


PARENT PROGRAM | 

Statistical Analysis of a Monte Carlo Transport 

02 14 

CATEGORY 

02 IS SITE 

02 18 PROGRAM NO. 


05 14 CONTACT 



Je J. 

LamMotte 

| OATES | 

05 50 INITIATED 


05 54 COMPLETED 

0967 


IO 69 


WHO TO CONTACT ABOUT THE PROGRAM 

105 28 SITE 1 05 31 ORGN CODE |CS 39 PROJECT NO 


05 28 SITE 05 31 ORGN CODE 

LAR 11.150 

Q5 58 REVISION CODE 

I \ A REVISION 
\ I B CANCELLATION 


KDP32l|- 


05 48 STATUS 

I 1 A. UNDER DEVELOPMENT 
1 1 B. OPERATIONAL 

C. COMPLETED 


05 49 

D a. THIS PROGRAM j 
IS NOT FOR 
SHARING 


TIME AND COST FOR DEVELOPMENT 

05~59 MANMONTHS 1 05^64 MACHINE 05 '69 COMPUTER TYPE 05 74 TOTAL COST 
HOURS (DOLLARS) 

4*0 2* 0 20 0 O * 0 

— 1 1 - 1 — 1 — 1 y. ^ n nn i i — i i — i — i 

59 60 61 62 63 64 65 66 67 68 OUUU 74 75 76 77 78 79 80 

ELITE MARGIN I PICA MARGIN 


This program was written in conjunction with the transpc 
carried out in R2092 , B2091, and E20Q5* It uses the pai 
Mst ori es on the tape generated from the transport to 
determine the particle density with respect to -position, 
MentMj a nd /or energy _at prescribecLdegths v ithi n.the 
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APPENDIX B 


PROGRAM FLOW CHARTS 


Separate flow charts are not given for both INTPOL1 and INTPOL2 since they would 
be identical. Similarly, any subprograms common to both TWOGEV and PROTOS which 
would have identical flow charts are presented only once as a part of the TWOGEV Flow 
Charts. 
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TWOGEV FLOW CHARTS 
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PROTOS FLOW CHARTS 
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SUPER B FLOW CHARTS 
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STAT FLOW CHARTS 
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TABLE I.- ELEMENTS AND MIXTURES ON MASTER DATA TAPE 


Program 

Element or mixture 

Symbol 

Atomic 

weight 

Atomic number, Z, 
for INTPOL 

Density 

No. atoms per 
cm3 x io“24 

TWOGEV 

Oxygen 


O 

16 

8 



Aluminum 

AL 

27 

13 




Carbon 


C 

12 

6 




Oxygen 


O 

16 

8 




Aluminum 

AL 

27 

13 




Chromium 

CR 

50 

24 



PROTOS 

Copper 


cu 

63 

29 




Ruthenium 

RU 

96 

44 




Cerium 


CE 

136 

58 




Tungsten 

W 

180 

74 




Lead 


PB 

204 

82 




Uranium 

U 

238 

92 




Tissue* 


c 27 H 21 ON 






Aluminum 

AL 

27 

13 




Carbon 


C 

12 

6 




Silicon 


SI 

28 

14 




Copper 


CU 

63 

29 




Lead 


PB 

204 

82 




Iron 


FE 

54 

26 




Water 

Hydrogen 




1.0 

0,067 



Oxygen 




.0335 


Polyethylene I 

Carbon 

CH2 



0.950 

0.04087 


Hydrogen 





.08174 



Carbon 





0.03936 


Polyethylene n 

Hydrogen 

CH2 



0,915 

.07873 

SUPER B 


C 





0.02083 



H 





.06683 


Nylon 

O 

C 12 H 22°2 N 2 



1.1399 

.006076 



N 





.06074 



C 



1 


0,04769 


Tissue* 

H 

O 

c 27 h 21 on 



1,0999 

,03710 

.00177 



N 





.00176 



C 





0.03581 


Lucite 

H 

c 5 h 8°2 



1.1892 

.05730 



O 





,01433 


Air 

N 




1.0 

0,0344 



O 




.007527 


"'Tissue is already in a form acceptable by PROTOS so that INTPOL2 will not be necessary. 
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Figure 6.- Upper and lower limits of dose rate in extremities from the prompt 
spectrum of the February 1956 solar event. 
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